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In this paper, we compute the accuracy at which the planned space-based gravitational wave 
detector LISA will be able to observe deviations from General Relativity. To do so, we introduce six 
correction parameters that account for modified gravity in the second post-Newtonian gravitational 
wave phase for inspiralling supermassive black hole binaries with spin precession on quasi-circular 
orbits. The precession of the spins and the angular momentum modulate the gravitational waveform, 
resulting in additional structure which could reduce correlations in the parameter space and increase 
the detection accuracy of the alternative theory parameters. Also, the use of higher harmonics could 
create further structure and increase the time during which the signal lasts in the frequency window 
of LISA. In order to find error distributions for the alternative theory parameters, we use the Fisher 
information formalism and carry out Monte Carlo simulations for 17 different binary black hole 
mass configurations in the range 1O 5 M0 < M < 1O S M0 with 10 3 randomly distributed points in 
the parameter space each, using the full (FWF) and restricted (RWF) version of the gravitational 
waveform. We find that the binaries can roughly be separated into two groups: one with low 
(< 10 7 M Q ) and one with high total masses (> 1O 7 M ). The RWF errors on the alternative theory 
parameters are two orders of magnitude higher than the FWF errors for high-mass binaries while 
almost comparable for low-mass binaries. Due to dilution of the available information, the accuracy 
of the binary parameters is reduced by factors of a few, except for the luminosity distance which is 
affected more seriously in the high-mass regime. As an application, we compute an optimal lower 
bound on the graviton mass which is increased by a factor of ~ 1.5 when using the FWF. 

PACS numbers: 04.30.Db, 04.50.Kd 



I. INTRODUCTION 

Although General Relativity (GR) has so far passed 
all experimental and observational tests [T], some unsat- 
isfactorily explained phenomena still remain which could 
be more elegantly described by alternative gravity theo- 
ries. Among these theories are the proposed inflationary 
epoch of the universe shortly after the big bang which 
explains the temperature homogeneity of the cosmic mi- 
crowave background, dark matter which should account 
for the missing 23% of the mass in the universe and 
dark energy introduced as an attempt to drive the ob- 
served late accelerated expansion of the universe. More- 
over, attempts to quantize GR or to unify gravitation 
with the other three fundamental forces have so far all 
failed. Consequently, several modifications to GR have 
been proposed. Certain alternative theories work by in- 
troducing additional fields to the Einstcin-Hilbert action 
of GR. Scalar-tensor field theories such as Brans-Dickc 
theory [5] are candidates for reproducing inflation. Mod- 
ified Newtonian Dynamics (MOND) 3J attempts to get 
rid of dark matter by modifying the 1 /r 2 behavior of the 
gravitational potential; a relativistic version introducing 
scalar and vector fields called Tensor- Vector- Scalar grav- 
ity (TeVeS) has also been proposed [4]. The class of f(R) 
theories [5] modify the Einstein-Hilbert action by replac- 
ing the Riemann scalar by a function of the Riemann 
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scalar. More phenomenological approaches such as Mas- 
sive Graviton theories [BJ [7] study the wave propagation 
of a massive gravitational field. Since alternative theories 
can be heavily constrained by the observation of Solar 
System effects and pulsar binaries [T], they are expected 
to reduce to GR in the limit of weak fields and to con- 
tribute significantly only in spacetime regions with strong 
dynamical gravity such as binary black holes (BBHs). A 
good review of currently discussed alternatives to GR can 
be found in the appendix of [8]. 

Among the most popular gravitational wave detectors 
are laser interferometers. Several ground-based interfer- 
ometers such as LIGO (USA), Virgo (Italy) and GEO600 
(Germany) have been built and are already operating, 
being sensitive to high frequencies between 10 Hz and 
1 kHz. Currently LIGO is being upgraded to Advanced 
LIGO with a sensitivity ten times better, and is expected 
to observe several events per year and make gravitational 
wave detection likely within the next five years. Hence 
gravitational waves could finally be observed directly a 
hundred years after their theoretical prediction by Ein- 
stein. 

Complementary to ground-based detectors restricted 
by their short arm-length and seismic noise at low fre- 
quencies, the spaceborne detector LISA (Laser Inter- 
ferometer Space Antenna) is planned, originally as an 
ESA/NASA collaboration. Scheduled to be launched in 
the early 2020 's and consisting of three spacecrafts sepa- 
rated by five million km in a triangular shape, LISA will 
form a giant laser interferometer sensitive to frequencies 
in the band between 10 -5 and 1 Hz. Recently, NASA 
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has discontinued their participation in the LISA project; 
the European Space Agency will try to realize the project 
on their own with a reduced, affordable budget [5]. The 
LISA science case is currently reviewed and several al- 
ternative candidates are proposed. In this paper we only 
discuss the originally planned LISA detector. 

Among the strongest sources which LISA will detect 
are supermassive black hole binaries with masses between 
10 5 — 10 7 A/ Q . After a long inspiral phase, such binaries 
could merge into one single Kerr black hole which rings 
down from its excited state by emitting gravitational ra- 
diation. Compact binary inspirals produce a very clean 
and long-lasting gravitational signal which may be accu- 
rately described by harmonics of the orbital phase using 
the post-Newtonian (PN) formalism. Inspiralling BBHs 
emit gravitational radiation carrying information about 
binary parameters such as the black hole masses and the 
individual spins in its amplitude and phase. By using 
matched filtering techniques |10[ lllj , the binary parame- 
ters can be extracted from the noisy signal measured by 
the detector. Alternative gravity theories will also leave 
their imprints on gravitational waves, since they modify 
the strong-field dynamics of the BBH, resulting in a dif- 
ferent orbital phase evolution. Also a possible graviton 
mass will influence gravitational waves on their way to 
us by making their velocity frequency dependent. Since 
alternative theories are heavily constrained and LISA is 
expected to observe signals with very high signal-to-noise 
ratio (SNR), a signal from a BBH will be detected with 
GR waveform templates regardless whether or not GR is 
true. This could create a fundamental bias [8] in param- 
eter extraction if the signal is fitted with an incorrect 
GR waveform template, leading to incorrect parameter 
estimation. To fix this bias, additional parameters con- 
trolling deviations from GR can be introduced. Adding 
parameters while having the same information from the 
detectors increases the correlation between the extracted 
parameters and thus decreases the accuracy in the recov- 
ered parameter values. 

Previous papers computed bounds which LISA could 
place on the Brans-Dicke parameter wbd (see e.g. [TH 
[T5] ) or on the graviton Compton wavelength X g (see e.g. 
[BJ) using matched filtering. Due to the no hair theo- 
rem, for BBHs, scalar field effects in Brans-Dicke theory 
arising from the inner structure of compact objects can- 
not be distinguished; however, such massive binaries are 
an excellent environment to test massive graviton effects. 
After a first analysis of massive graviton propagation by 
Will [5J, Berti et al. Q3] introduced spin parameters 
and spin-orbit/spin-spin couplings, finding a loss of accu- 
racy due to the extra parameters included in the model. 
Stavridis and Will [7] considered the full precession of 
the spins and discovered that the resulting phase modu- 
lation restores the lost accuracy on A g . Yagi and Tanaka 
[15j included eccentricity to the system and found that 
the additional structure through both precession and ec- 
centricity increases the measurement accuracy by an or- 
der of magnitude. Arun and Will [TB] showed that the 



bounds on X g are improved by almost an order of mag- 
nitude for non-spinning BBHs when using the full wave- 
form (FWF) instead of the restricted waveform (RWF) 
which takes the phase up to full PN order but consid- 
ers the amplitude only to leading order. Taking higher 
harmonics into consideration increases the time during 
which the signal stays in the frequency window of LISA 
and shows a richer structure in the gravitational wave, 
leading to less correlation in the parameter space. Kep- 
pel and Ajith |17j used hybrid inspiral-merger-ringdown 
waveforms and found that they lead to a ~ 10 times 
higher accuracy than for inspiral-only waveforms. More- 
over, Berti et al. [H] pointed out that the combination 
of the bounds on X g from individually observed inspi- 
rals in a two-year running time can again raise the ac- 
curacy by an order of magnitude. Tables summarizing 
lower bounds on X g and upper bounds on wbd found 
by previous works are e.g. provided by |15[ I17j. Arun 
et al. PTSJ re-interpreted the matched filtering method 
and fitted the post-Newtonian coefficients to the wave- 
form instead of the parameters usually extracted from 
them. They discussed to what extent LISA will be able 
to measure deviations from the 3.5PN gravitational wave 
phase parameters in General Relativity. Yunes and Pre- 
torius [8] generalized this approach to a parameterized 
post-Einsteinian (ppE) framework which maps different 
types of alternative theories to the gravitational wave- 
form of a compact binary merger. Cornish et al. |20j 
used Markov Chain Monte Carlo simulations to investi- 
gate parameter biases and possible bounds on the ppE 
parameters. 

In this work we parameterize alternative theories by in- 
troducing corrections to the post-Newtonian coefficients 
of the orbital phase for a BBH inspiral, including the 
full 2PN precession of spins and angular momentum. We 
use the FWF up to 2PN instead of the RWF, summing 
first six harmonics. We postpone the discussion of eccen- 
tric orbits to later work and restrict our calculations to 
quasi-circular orbits. Since matched filtering is far more 
sensitive to the gravitational wave phase than to the am- 
plitude, we do not consider corrections to the amplitude 
of the wave. We evaluate the measurement accuracy with 
which standard LISA will be able to detect such correc- 
tions for BBHs. To estimate the errors on the param- 
eters, we make use of the Fisher information formalism 
which is legitimate in the limit of high SNR which LISA 
will provide. 

The organization of this paper is as follows. In Sec. 
im we shortly introduce the necessary equations to de- 
scribe the evolution of the inspiral phase, the spins and 



the angular momentum of a BBH up to 2PN. In Sec. Ill 
we introduce small departures from GR into the post- 
Newtonian frequency evolution equation. We then com- 
pute the modified orbital phase evolution in this scheme, 
incorporate it into a modified waveform template in Sec. 
WJ taking the waveform to be the sum of harmonics of 
the orbital phase, and compute the Fourier transformed 
waveform including alternative theory parameters. In 



3 



Sec. |V] we review the Fisher information formalism in 
order to estimate the errors on the parameters. In sec. 
| VI| we explain the details of the Monte Carlo simulations 
we carried out. We discuss the resulting error distribu- 
tions on selected parameters in Sec. 



VII to see to what 



extent we can measure deviations from the 2PN gravita- 
tional wave phase predicted by GR and how strongly the 
binary parameters are affected by the introduction of six 
new parameters to the model. We discuss two represen- 
tative BBH systems in Sees. VII A I and VII B I. Because 



systems at higher redshifts experience higher errors, we 
plot the maximal redshifts for different upper error lim- 



its of the alternative theory parameters in Section VII C 



As an example we calculate the resulting optimal lower 
bounds on the Compton wavelength of the graviton in 
Sec. ~ 



VII D 



We summarize our work and discuss possi- 
ble extensions in Sec. VIII In Appendix [A| we give the 
expressions we used for the 2.5PN and 3PN frequency 
evolution. We give tables with best-case, worst-case and 
median measurement errors of both the binary and alter- 
native theory parameters in Appendix |B| 
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is a dimensionless orbital frequency parameter, M — 
mi + m2 is the total mass and v = vn,\vn,2 /M is the 
symmetric mass ratio. The spin-orbit and spin-spin cou- 
plings are given by 
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II. EVOLUTION OF BLACK HOLE BINARIES 
WITH PRECESSING SPINS 
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A complete description of the inspiral evolution of two 
spinning black holes on a quasi-circular orbit with two in- 
dividual masses mi 2 an d the corresponding spin vectors 
Si, 2 (t) is given by the angular momentum unit vector 
L(t), the orbital angular frequency u)(t) and initial val- 
ues for the orbital phase tp(to) and frequency u)(to). Fur- 
ther characteristics such as the orbital separation can be 
related to oj using post-Newtonian expressions. There- 
fore a quasi-circular BBH inspiral can be described by 12 
intrinsic parameters. In order to relate the binary with 
a detector, a unit vector h pointing from the detector 
to the barycenter, and a luminosity distance d^ between 
the two can be introduced, bringing an additional set 
of 3 extrinsic parameters into play. Thus, to describe a 
BBH inspiral on quasi-circular orbit, 15 parameters are 
required. 

Since a description of the motion of such a system with 
full General Relativity is only possible with numerical 
methods and at high computational cost, an analytic ex- 
pansion of the Einstein equations in powers of v/c has 
been studied: the post-Newtonian (PN) formalism. Cur- 
rently, the equations of motion for spinning objects are 
known up to 2.5PN, while spin-spin and spin-orbit cou- 
pling terms are only known up to 2PN |21) . Therefore we 
take all the relevant expressions up to 2PN, i.e. 0[(v/c) 4 ] 
away from leading order. The evolution equation for the 
angular frequency of a BBH system is [2"2"] 



respectively. The precession of L and S\ t 2 induces a 
time dependence for these couplings, and thus a modula- 
tion of the gravitational wave phase. Such a modulation 
is good for parameter estimation, since it creates extra 
structure which helps us decorrelate the entangled pa- 
rameters. The orbit-averaged evolution equations with- 
out radiation reaction (L + Si + S2 = 0) at leading PN 
order are [231 



a_ 

c 2 r 



2 + 



2m 2 



x L 



3G 
2c2 



Si = 



$ [( 



^((S 2 L|S, + (S, . • L)S.) L 

3m,- \ , 1 c , 3 



2 + 



2 v Si • L 



(5) 

) £] x Si, (6) 



with i 7^ j and i,j € {1,2}. The orbital separation r 
and the angular momentum are related to the orbital 
frequency by the Newtonian relations 
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since higher-order corrections would exceed the 2PN or- 
der. Eqs. ( 1 ) and Q enable us to express the evolution 
equations (6) in terms of the frequency u>: 
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We express the gravitational wave phase in terms of 
the "principal + direction" [24 defined as the direction 
of the vector L x n. A precession of the angular mo- 
mentum vector changes the principal + direction. The 
resulting modulation of the gravitational waveform can 
be expressed by modifying the phase by 
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where ojq is the orbital frequency at time to, 5ipo 



— j t ° (dS(p/ dt)dt, and dL/dui is given in Eq. (10). The 
resulting 2PN orbital phase is then, expressed in terms 
of the orbital angular frequency: 4>(uj) — (f{uS) + 5tp(u>). 

A signal observed from a BBH at cosmological distance 
is redshifted, i.e. the observed frequency is f a = / c /(l + 
z), where f e is the frequency of the gravitational waves 
emitted by the binary. The relation between redshift 
and luminosity distance in a ACDM cosmology without 
radiation and with f^A = 0.72, tt m = 0.28 and H = 70.1 
km/s/Mpc [23] is 



, . . . x c f dz' 
d L (z) = (l + z) w =■ (14) 

The previous derivations have been done for a BBH 
at z = 0. For binaries at cosmological distance, the 
masses are redshifted according to mi, 2 = (1 + z)mi^ 
and the distance between the observed system and the 
observer must be replaced by the luminosity distance. 



Unfortunately for gravitational wave experiments, it is 
not possible to disentangle redshift, mass and distance: 
only two parameters out of these three can be inferred. 
Simultaneous observations of electromagnetic counter- 
parts, through which the actual redshift could be mea- 
sured, could break this correlation and lead to interesting 
astrophysical insights. 



III. MODIFICATIONS TO THE 2PN ORBITAL 
PHASE 

Matched filtering techniques are more sensitive to the 
gravitational wave phase than to the amplitude. The 
signal from a BBH inspiral can be described as a sum 
of harmonics of its orbital phase; to find the imprints of 
alternative gravity theories on gravitational waves it is 
therefore reasonable to look at how the orbital phase evo- 
lution of a BBH changes for small departures from GR. 
In the 2PN expansion, the orbital phase evolution can 
be found by integrating the frequency evolution equation 
(see Eq. for the PN coefficients bi) 
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As thoroughly discussed by Yunes and Pretorius in the 
derivation of their ppE framework [S], in the adiabatic 
approximation the dimensionless frequency can be ex- 
pressed as 
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E is the total binding energy or Hamiltonian (conser- 
vative part) of the system while E stands for the en- 
ergy loss through gravitational waves or other physical 
degrees of freedom of energy loss (dissipative part) . Con- 
sidering the impact of alternative theories on these two 
quantities leads to modifications of the gravitational wave 
phase. Certain theories such as Brans-Dicke theory in- 
troduce scalar fields which lead to a difference in the 
self-gravitational binding energy Q per unit mass [26] . 
producing additional dipole radiation. The energy loss 
formula including dipole contributions can be expressed 
to leading quadrupole order as [5] 157] : 
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Here, v and r are the orbital velocity and separation of 
the system, respectively, while n\ and k 2 are so-called 
Peter-Mathews parameters and is a coefficient for 
the dipole contribution. £ thcr stands for any other 
energy loss channel, either through other polarizations 
or as yet unknown physical processes. Since we do not 
have any good parametrization for £ Q thor so far, we do 
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not consider it. In terms of dimensionless frequency, the 



dipole radiation term in Eq. (17) leads to an additional 



term in the PN expansion ( 15 ) 



We introduce a general paramctrization where the 
effects on the phase are emphasized and no corrections 
to the wave amplitude are considered. The calculations 
are done for quasi-circular binaries with precession 
of both black hole spins described by the full 2PN 
waveform (2PN expansion of both the phase and the 
amplitude). We start by introducing corrections to the 
2PN orbital frequency evolution dx/dt which will lead 
to a corrected version of the 2PN orbital phase. To do 
that, we introduce a correction term a,i for every 2PN 
coefficient bi and an additional x~ x and x 1 / 2 term. The 
waveform is a series of harmonics of the orbital fre- 
quency: this results in corrections to both the amplitude 
and phase of the gravitational wave in frequency space. 
Its derivatives give the Fisher matrix elements needed 
for parameter estimation. Products of a correction 
term a^ix^ 1 with a PN expanded expression such as 
l + bxx + b. 



'3/2' 



^3/2 _|_ resu i|; m b 2 featuring already at 
1PN order. Hence for the final result to be consistent at 
2PN order, we need to do all the calculations up to 3PN, 
truncating at 2PN only at the very end. The current 
2.5PN expansion accounts for spin-orbit effects while 
the 3PN expansion does not consider spin effects at all. 
Nevertheless, these higher order expansions can be used 
as approximations. The 3PN evolution equations of the 
dimensionless orbital angular frequency are, motivated 
from H21HB] (see appendix \M 



3PN 



64i> c 3 



1 + b x x + b 3/2 x 3/2 + b 2 x 2 



+b 5/2 x 5/2 + b 3 x 3 + 6 3 ,io g a; 3 logO) 



(18) 



h = 

h/2 = 

b 2 = 

h/2 = 

b 3 = 

\ log = 



/743 llv 
V 336 + ~4~ 



4^-^8(113,75) 



/34103 13661^ 59^ 2 



V 18144 2016 
4159 189f 



672 

x/3(l,0) + 



18 
1 

c 



- — ct(247, 721' 
48 y 

40127 1465^ 



1008 



28 



16447322263 
139708800 



583 
1712 7e 



3049i/ 
168 
16tt 2 



£(-M) , 
56198689^ 



105 



217728 



451tt 2 i/ 541i/ 2 



48 
856 

105' 



896 



5605i/ d 
2592 



856 
105 



log(16), 



(19) 



where (3 and a are the spin-orbit and spin-spin couplings, 
respectively. To account for alternative theories we gen- 
eralize the frequency evolution to 
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including corrections to every existing PN parameter 
and an additional x~ x and x 1 ^ 2 term. The reason why 
x 2 log(a;) appears is that a term proportional to x 3 log(x) 
enters the 3PN phase which has to be included in 2PN 
corrections because of couplings with x~ x terms. 

We now follow the steps for the derivation of the grav- 
itational waveform presented in |29j . introducing these 
additional corrections, keeping them at first order, and 
truncating at 3PN. 



By inverting and integrating Eq. ( 20 1 we find the time 



t(x) as a function of the frequency to be of the form: 
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with 



To find the orbital phase as function of frequency, we 
need to recast t(x) into a series expansion for x(t); we 
are then able to find the phase by integrating uj oc x 3 ! 2 
over time: 
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with the phase corrections Ai as functions of the original 
orbital frequency evolution corrections a^. At this point 
we choose not to consider the correction term A 2 ^ og in 
our framework for simplicity and thus set A 2 j og — in 
the following. 



IV. MODIFICATIONS TO THE 2PN 
WAVEFORM 

Having found a 2PN expression for the orbital phase 
corrections, we are able to construct the gravitational 
waveform as a series of harmonics of the orbital fre- 
quency: 
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Here, (j> is the orbital phase of the binary with spin preces- 
sion included: cj)(t) = [<f{t)] mod + Sip(t). The coefficients 

-A+x' ^+"x are both post-Newtonian series in x: 



On and <Pn are the spherical angles of the position of the 
binary in the detector frame, and 4>n is defined through 
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with ip = [<p] mod + Sep + (f>D, including the LISA Doppler 
phase (f>D(t) = (wR/c) siuOn cos($(t) — 4>n), where R — 
1 AU and 4>(t) = 2irb/\ yr as explained in [29 . The 
harmonic coefficients are 
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By changing the cosine+sine representation into a co- 
sine+phase representation, we can write Eq. (25 1 as 
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Explicit expressions for A+ x and can be found 

in [29]. LISA will form two different detectors with un- 
corrected noise: for a detector k with antenna pattern 
functions and Fj* , the response function is 
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The Fourier transform of the response function is then, 
writing the cosine as an exponential and defining the new 
phase tp kin = n([(p] mod + 5tp + cj> D ) + <j) p k ol n . 
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The n = integral accumulates around frequencies 
different from the gravitational wave frequency and 
e i(27r/t+i> fc- „) arounc j negative frequencies, so both can be 
neglected. Then the Fourier transform reduces to 
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In the stationary phase approximation (SPA, see e.g. 
I3"T]L hk{f) is approximated by 
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evaluated at the stationary points t n = t 2 p^s(f /n). The 
square root of the reciprocal of the second derivative of 
ipk,n is found to be 
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+S 2 x 2 + S 2i i og x 2 log(a;). 



(38) 



The Si are functions of the orbital phase corrections Ai. 
The waveform can then be written as 
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+ *3/2^ 3/2 + ^2X 2 



(40) 



The are also functions of the orbital phase corrections 
Ai. It makes thus sense to work only with the phase 
correction parameters ^i from now on. The correction 
coefficients of AS are then, given as functions of Vf^: 
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M/) ~ f E A^(t(f/n))x-^[S(f/n)] mod 

n>l 

— 6ip(f/n) 
-Mt(f/n)})-^[t(f/n)}}}, 

(39) 

where the modified phase is defined as [*(/ /n)] mod = 
[*(//n)]apN + A* , with 



V. PARAMETER ESTIMATION 

To estimate how accurately LISA can measure devia- 
tions from the 2PN gravitational wave phase predicted 
by General Relativity, we use the standard Fisher in- 
formation formalism for gravitational wave experiments, 
as reviewed in [32l [33] . The Fisher information formal- 
ism holds only in the limit of high SNR; this is true for 
LISA for which SNR of a few thousands are expected. 
For low SNR, advanced Bayesian techniques exploring 



the whole parameter space such as Markov Chain Monte 
Carlo methods, (see e.g. [2Ql[34]) are needed. Also, once 
data will become available, Bayesian statistics taking into 
account prior probability distributions will be the pre- 
ferred framework. 

We assume the gravitational wave signal to be buried 
in stationary Gaussian noise n(t) such that the different 
Fourier components n(f) are uncorrelated. Moreover, we 
presume that the noise of the two detectors is totally un- 
correlated. For a signal h(t) described by a true param- 
eter set t , with noise with spectral density S n (f), the 
probability for the measured data d(t) = n(t) + h(t;6 t ) 
to take this specific form is proportional to 
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), (49) 



where L = 5 x 10 9 m is the arm length of standard 
LISA, S p = 4 x 10 -22 m 2 Hz -1 is the white posi- 
tion noise level, S a = 9 x 10 -30 m 2 s~ 4 Hz -1 is the 
white acceleration noise level, and /* = c/(2ttL) is the 
arm transfer frequency. The total noise curve is then 
S h (f) = S n (f) + S conf (/). 



p{d\6 t ) oc e -{d-h{e t )\d-h{e t) )^ 
where the inner product (g\h) is defined as 



(42) 



w= 4 ^r r Sr d/ - (43) 

The use of a waveform template with the parameter set 
6 is inaccurate by A6 l = B\ — Q % . The errors A9 1 are then 
approximately given by maximizing the above likelihood 
distribution, expanding it in the errors assumed to be 
small and keeping only first derivatives [33] : 



where S is the covariance matrix and 



(44) 



dh 


dh 


d¥ 





(45) 



is the so-called Fisher matrix. The expected measure- 
ment errors on the parameters 9 l can be expressed as 



A6 l = ^(r- 1 )". 



(46) 



We chose the same noise curve for standard LISA as in 
[2"9"] . namely the piecewise fit used by the LISA parame- 
ter estimation community [35) given by the instrumental 
noise 
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VI. SIMULATIONS 

For our simulations 21 parameters are needed: 15 GR 
parameters plus 6 alternative theory parameters. We use 

(i) log lo mi/M and log 10 m 2 /M Q , for the masses of 
the two black holes. 

(ii) Hi = cos 9i and <f>i, for the spherical angles of the 
orbital angular momentum L at 7 = g. 

(iii) Hi = cos 9\ and <pi for the spherical angles of the 
spin of the first black hole Si at 7 = | . 



(iv) X 



\Si\ for the dimensionless strength of 



1 — Gra\ 

the spin of the first black hole, which has to satisfy 
0<xi<l- 

(v) /12 = cos 62, 4>2, and X2 for the second black hole, 
defined equivalently as for the first one. 

(vi) t c , the time of coalescence. 

(vii) (p c , the phase at coalescence. As this phase is ran- 
dom and its determination is not of any astrophys- 
ical interest, we can safely neglect constants in the 
orbital phase, in particular Sipo from Eq. (13 1. 



(viii) ji n = cos#„ and <j) n , the spherical angles of the 
position of the binary in the sky. 

(ix) e?£, the luminosity distance between the source and 
the Solar System. 

(x) * 4 with i e {-1, 0, 1/2, 1, 3/2, 2}, the 6 alternative 



theory parameters defined in section IV 



and the confusion noise 



All angles are taken in the frame tied to the distant 
stars. Moreover, we set t = to be at the time when 
LISA will start operating. 

We perform Monte Carlo simulations, keeping the 
masses the redshift z and the alternative theory pa- 
rameters fixed, and randomizing all other parameters 
using a flat probability distribution. The spin precession 
equations Q are integrated using a fourth order adap- 
tive Runge-Kutta algorithm to find the evolution of L(uS) 
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and Si t 2(u}), going backwards in frequency. As a good 
starting point for u we chose the frequency at the inner- 
most circular orbit (ISCO) for a Schwarzschild black hole: 
W2pn (nsco = 3 r s) wnere r s is the Schwarzschild radius, 
as pointed out in [55]. We stop the evolution either at 
t = or when the frequency of the highest harmonic goes 
below the LISA band (6w < 3 x 1CT 5 Hz). The upper 
and lower bounds on all the randomized parameters of 
the simulation are straightforward (cLl is just a function 
of the redshift z, defined in (14)), except for t c for which 
we set a lower bound of t c = t 2 -p-^(ui{r — rigco)) us- 
ing Eq. (21) and an upper bound of t c = 2yr, which is 



the minimum science requirement for the LISA mission 
running time. 

Using the angular momentum, spin and orbital time 
evolution we are able to compute the Fisher matrix ( 45 ) 
taking the analytical derivatives with respect to t c , d,L, 
4> c , n n , <p n and all the GR correction parameters The 
first three derivatives are easy to compute: 



d~h k (e j ,f) 

dtr. 



dh k (e\f) h k (&,f) 



dd L 

dh k (0 3 J) 
dip c 



'^nhk, n (0 J ,f), 



(50) 
(51) 
(52) 



where hk, n is the nth harmonic of h k . The derivatives 
with respect to the corrections are of the form 



avoid matrix inversion problems, we use a normalization 
of the Fisher-Matrix so that all diagonal elements are 
An — 1 and all off-diagonal elements are in the range 
Aij €[-!;!]: 



(55) 



v/r~r 



j j 



After inversion, the covariance matrix can then be recov- 
ered with 



AT 1 



(56) 



OJ 



In situations where L ■ n is close to 1, the Runge Kutta 
method fails to converge because 



dS(p 



1- lL-n 



Whenever this happens, we take the approximate value 



6<p(<jJ + Sid) R: 

Sip(u)) + angle 
as explained in 
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(53) 



and can be calculated in a straightforward way. The 
derivatives which we could not compute analytically are 
approximated by 



dhkj^J) ^ h k (0 j + fg_A /) - h k (0i - e^/2, /) 
d6 l e 

(54) 

where e is a small displacement of the parameter 0j which 
we chose to be of the constant value e = 1CP 7 for every 
parameter, except for (pi for which e was divided by 2 — 
2|w|, A*i (i € {1, 2}) for which e was divided by 5\i, and 
(pi for which e was divided by 10xi(l — \fH\)- The formula 
is accurate up to 0(e 2 ). 

For each set of parameters we then compute the Fisher 
matrix using Clenshaw- Curtis quadrature and then in- 
vert it in order to find the corresponding errors on the 



parameters which we analyze in section VII In order to 
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VII. RESULTS 



RWF21 
RWF15 




10"* 



0.001 



0.01 



0.1 



Am, / m. 



FIG. 1: Comparison between the estimated distributions of 
the measurement error on mi for a low-mass binary system 
mi = 1 X 10 6 M© and m 2 = 3 x 1O 5 M with (RWF21) and 
without (RWF15) including alternative theory parameters, 
using only the restricted waveform. 

We performed simulations for 17 different mass config- 
urations, with total masses between 10 5 M Q and 10 s M , 
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mass ratios varying between 1:1 and 1:10, and using 10 3 
points in the parameter space for each configuration. 

The redshift has been kept fixed to z = 1 since it is 
not possible to disentangle redshift, mass and distance. 
The signal coming from a binary with masses mi. 2 at 
redshift z and luminosity distance cIl(z) can be expressed 
with one from an apparent binary with mi. 2 = T+jj^ m i,2 
at redshift zq and luminosity distance cIl(zo) multiplied 
by an overall factor of dz(zo)/d,L(z). Thus every BBH 
inspiral producing a signal at redshift z can be described 
with a waveform template at redshift zq. The Fisher 
matrix scales as 



r (z) = 

ij 



dh 



(m 1 ,m 2 , z) 



dh 



(mi,m 2 ,z" 



<Il(zo) 
d L (z) 

dL(zp) 
dL{z) 



dh 

^-(mi,m 2 ,zo) 



-<( z o) 



dh 



(mi,m 2 , z ' 



The errors on the parameters scale then with 



A6\z) = ^AFiza). 
dL(zo) 



(59) 



(60) 



Since we choose to work in a picture where General 
Relativity is the theory assumed to be true and we are 
keen to know how well LISA will be able to measure 
deviations from its post-Newtonian expansion terms ?/>i, 
we fixed the alternative theory parameters to the fiducial 
values 9i = 0. 

For each of the 17 binaries we computed the best-case 
measurement error (5% quantile), the typical error (me- 
dian) and the worst-case error (95% quantile) for the 



full (FWF) and restricted waveforms (RWF) and present 
them in tables l] |XIV For each BBH parameter we are 
interested in, we give an error table with (21 parameters 
in total) and without (15 parameters in total) includ- 
ing the alternative theory parameters We do this 
to show how much accuracy is lost by introducing alter- 
native theory corrections into a GR waveform template. 
For binaries where no signal can be extracted from the 
dataset, we fix the error to infinity. 

We give the errors on the sky localization not in terms 
of errors on \i n and <p n but instead in terms of an error 
ellipse with principal axes 2a and 2b, enclosing the region 
outside of which there is an l/e probability of finding the 
binary, following [36] . 

Moreover, in tables |XVj|XX| we give measurement er- 
rors on the alternative theory parameters, using both the 
RWF and FWF. 

We roughly divide the binaries into two classes: low- 
mass binaries (M < 10 7 M Q ) and high -mass binaries 
(M > 1O 7 M ). Below we discuss these two cases, us- 
ing BBHs with mi = 1O 6 M , m 2 = 3 x 10 5 M Q and 
mi = 3 x 1O 7 M0, m,2 = 10 7 Af Q as representative ex- 
amples for low-mass and high-mass binaries respectively. 
We find when using both the RWF and the FWF, the 
error distributions of the mass and spin parameters be- 
have similarly, losing a factor 1.2 — 5 of accuracy when 
alternative theory parameters are included. The error on 
the orientation of the binary 2a and 2b is at maximum 
an order of magnitude worse. For high-mass binaries, 
factors of ~ 10 and ~ 100 are lost in the determination 
of the luminosity distance d^, using the FWF and RWF 
respectively. While the RWF/FWF errors on the alter- 
native theory parameters are almost equal for low-mass 
binaries, the RWF errors are about 100 times higher for 
high-mass binaries 
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FIG. 2: Comparison between the estimated distributions of 
the measurement error on mi for a low-mass binary system 
mi = 1 x 1O 6 M and m 2 = 3 x 10 5 A'/ Q with (FWF21) 
and without (FWF15) including alternative theory param- 
eters and using the full waveform. 



For low-mass binaries with total masses below 1O 7 M 
we find that in general, using the FWF instead of the 
RWF improves the measurement errors A\&j on the al- 
ternative theory parameters by a factor of ~ 1.5 — 3. The 
correlation with the new parameters causes a decrease in 
the accuracy of the 15 existing parameters. For both 
the FWF and the RWF, the errors on the mass and spin 
parameters mi, xi an d X2 are typically worse by a 
factor of 2 — 5 while the luminosity distance dz, is approx- 
imately half as accurate. The angular orientation errors 
2a and 2b increase only by ~ 10%; this is reasonable, 
since we do not expect alternative theories to correlate 
strongly with rotations on a large scale. Therefore it is 
not necessary to use the FWF instead of the RWF for the 
sole purpose of measuring alternative gravity parameters 
in the low-mass regime. 

We present selected distributions of the measurement 
errors Ami /mi, Axi/xi; 2a, Ad^/d^ and all the six 
A^i in figures [TpH The error distributions of Am 2 /m 2 , 
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^■Xilxi and 2b are similar to the ones of Ami /mi, 
Axi/Xi and 2a. 




FIG. 3: Comparison between the estimated distributions of 
the measurement error on xi for a low-mass binary system 
mi = 1 x 1O 6 M and m 2 = 3 x 1O 5 M with (RWF21) and 
without (RWF15) including alternative theory parameters, 
using only the restricted waveform. 
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FIG. 4: Comparison between the estimated distributions of 
the measurement error on xi for a low-mass binary system 
mi = 1 x 1O 6 M and m 2 = 3 x 10 5 M Q with (FWF21) 
and without (FWF15) including alternative theory param- 
eters and using the full waveform. 



B. High-mass binaries 

By using the FWF instead of the RWF for high-mass 
binaries with total masses > 10 7 M Q , we find significant 
improvements for the measurement errors of the alter- 
native theory parameters by factors of <~ 100 — 1000 for 
A*_i, - 30-60 for A* and A* 1/2 , and - 10-100 for 
A'J'i, A5 , 3/2 and A^ 2 . This makes it clear that it is in- 
evitable to use the FWF in the high-mass regime. In any 
case, since the second harmonic spends only a few orbits 
in the LISA band, the use of the RWF is not trustwor- 
thy. Moreover, for BBHs with total masses higher than 
10 8 M Q , LISA will not be able to see the second harmonic 




FIG. 5: Comparison between the estimated distributions of 
the major axis of the positioning error ellipse for a low-mass 
binary system mi = 1 x 1O 6 M and m 2 = 3 x 10 5 M Q with 
(RWF21) and without (RWF15) including alternative theory 
parameters, using only the restricted waveform. 
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FIG. 6: Comparison between the estimated distributions of 
the major axis of the positioning error ellipse for a low-mass 
binary system mi = 1 x 1O 6 M0 and m 2 = 3 x 10 5 Mq with 
(FWF21) and without (FWF15) including alternative theory 
parameters and using the full waveform. 
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FIG. 7: Comparison between the estimated distributions of 
the measurement error on (Ll for a low-mass binary system 
mi = 1 x 10 6 M Q and m 2 = 3 x 1O 5 M with (RWF21) and 
without (RWF15) including alternative theory parameters, 
using only the restricted waveform. 
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FIG. 8: Comparison between the estimated distributions of 
the measurement error on dh for a low-mass binary system 
mi = 1 x 10 6 M Q and m 2 = 3 x 10 5 M Q with (FWF21) and 
without (FWF15) including alternative theory parameters, 
using only the restricted waveform. 



FIG. 11: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
\]>i/2 for a low-mass binary system mi = 1 x 1O 6 M and 
m 2 = 3 x 1O 5 M , using the RWF and the FWF. 
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FIG. 9: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
for a low-mass binary system mi = 1 x 1O 6 M and 
m 2 = 3 x 10 5 M Q , using the restricted waveform (RWF) and 
the full waveform (FWF). 



FIG. 12: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
$i for a low-mass binary system mi = 1 x 10 6 M© and m 2 = 
3 x 10 5 M Q , using the RWF and the FWF. 
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FIG. 10: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
v&o for a low-mass binary system mi = 1 x 10 M© and m 2 = 
3 x 10 5 M Q , using the RWF and the FWF. 



FIG. 13: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
*3/ 2 for a low-mass binary system mi = 1 x 10 6 M© and 
m 2 = 3 x 10 5 M Q , using the RWF and the FWF. 
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FIG. 14: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
^2 for a low-mass binary system mi = 1 x 1O 6 M0 and m-2 = 
3 x 10 5 M Q , using the RWF and the FWF. 
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FIG. 17: Comparison between the estimated distributions of 
the major axis of the positioning error ellipse for a high-mass 
binary system mi = 3 x 10 7 M Q and m2 = 1 X 10 7 A-/q with 
(FWF21) and without (FWF15) including alternative theory 
parameters and using the full waveform. 
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FIG. 15: Comparison between the estimated distributions of 
the measurement error on mi for a high-mass binary sys- 
tem mi = 3 x 10 7 M Q and m 2 = 1 x 10 7 M o with (FWF21) 
and without (FWF15) including alternative theory parame- 
ters and using the full waveform. 
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FIG. 16: Comparison between the estimated distributions 
of the measurement error on xi f° r a high-mass binary sys- 
tem mi = 3 x 10 7 M Q and m 2 = 1 x 10 7 M o with (FWF21) 
and without (FWF15) including alternative theory parame- 
ters and using the full waveform. 



at all and so the RWF cannot be used. For both the FWF 
and the RWF, the errors on the mass and spin parame- 
ters mi, m2, Xl an d X2 are typically worse by a factor 
of 1.2 — 4 when accounting for alternative gravity pa- 
rameters. The luminosity distance is about 50 — 1000 
times less accurate for the RWF while for the FWF it is 
only ~ 10 — 100 times worse. For the FWF, the angular 
orientation error 2a and 2b is at maximum 5 times worse 
while the RWF loses up to a factor of ~ 10 in accuracy. 

We present selected distributions of the measurement 
errors Am^mj, ^Xi/Xii 2a, Ad^/d^ and all the six 
in figures [15][24] 
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FIG. 18: Comparison between the estimated distributions of 
the measurement error on 4 for a high-mass binary system 
mi = 3 x 10 7 M Q and m 2 = 1 x 10 7 Af Q with (FWF21) and 
without (FWF15) including alternative theory parameters, 
using only the restricted waveform. 
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FIG. 19: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
for a high-mass binary system mi = 3 x 1O 7 M and 
m 2 = 1 x 1O 7 M , using the RWF and the FWF. 



FIG. 22: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
for a high-mass binary system mi = 3 x 10 7 M Q and m 2 = 
1 x 10 7 M Q , using the RWF and the FWF. 
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FIG. 20: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
*o for a high-mass binary system mi = 3 x 10 7 M Q and m 2 = 
1 x 10 7 M Q , using the RWF and the FWF. 



FIG. 23: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
^3/2 for a high-mass binary system mi = 3 x 1O 7 M0 and 
m 2 = 1 x 10 7 M Q , using the RWF and the FWF. 
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FIG. 21: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
\Pi/2 for a high-mass binary system mi = 3 x 1O 7 M and 
m 2 = 1 x 10 7 M Q , using the RWF and the FWF. 



FIG. 24: Comparison between the estimated distributions of 
the measurement error on the alternative theory parameter 
<E' 2 for a high-mass binary system mi = 3 x 10 7 M Q and m 2 = 
1 x 10 7 M Q , using the RWF and the FWF. 
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C. Upper limits for redshifted masses 

All the errors tabularized in appendix [B] are given for 
the fixed redshift z — 1. Some of them in the high-mass 
regime are apparently too high at z = 1. Nevertheless, 
since the measurement accuracy of the parameters is cor- 
related with the redshift as given in Eq. (60 1, for an 



equivalent mass configuration at a lower redshift the er- 
rors should reduce to reasonable values. Since we do not 
know the actual values of the alternative theory param- 
eters, we cannot fix the accuracy with which we want to 
measure ^i(z). For this reason, we introduce the relative 
accuracy parameter a such that A^P^z)/^ < a where 
ipi is the fiducial 2PN phase coefficient from ^pn in Eq. 



( 40 1 . The maximal redshift is then given as 



(ad L {z ) |A*i(*b)M 



(61) 



where z(c?l) is the inverse of ( |14[ ) and can be com- 
puted numerically. We use here the 5%-quantilc for 
A\l/j(zo = 1) as given in the tables in appendix [B] i.e. 
we define the (optimistic) maximal redshift as the red- 
shift where 5% of the binaries in the sample can still be 
seen with relative accuracy less than a. Since we expect 
corrections to the 2PN phase parameters of GR to be 
small (at least for the lower PN orders), we focus here 
on a relative accuracy below 10%. At redshift z = 1 this 
accuracy is already difficult to reach for binaries with 
masses above 10 6 M Q (see also [19]). It is important to 
emphasize that we concentrate here on actually measur- 
ing the alternative theory parameters instead of just set- 
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FIG. 25: Maximal redshifts for the alternative theory param- 
eter such that the relative error A^o/V'o is smaller than 
a. ipo is the corresponding fiducial 2PN phase coefficient. For 
a relative error of 1%, low-mass binaries are suitable up to 
redshifts z ~ 1 — 10 while high-mass binaries can be observed 
up to z ~ 0.01 - 0.1. 
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FIG. 27: Maximal redshifts for the alternative theory param- 
eter ^3/2 such that the relative error A^ 3 / 2 /V , 3/2 is smaller 
than a. ipi is the corresponding fiducial 2PN phase coefficient. 
For a relative error of 1%, low-mass binaries are suitable up to 
redshifts z ~ 0.1 — 1 while high-mass binaries can be observed 
up to z ~ 0.01 - 0.1. 
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FIG. 26: Maximal redshifts for the alternative theory param- 
eter such that the relative error A^i/i/>i is smaller than a. 

is the corresponding fiducial 2PN phase coefficient. For a 
relative error of 1%, low-mass binaries are suitable up to red- 
shifts z ~ 0.1 — 1 while high-mass binaries can be observed 
up to z ~ 0.01 - 0.1. 



FIG. 28: Maximal redshifts for the alternative theory param- 
eter ^2 such that the relative error A^/V^ is smaller than a. 
ip2 is the corresponding fiducial 2PN phase coefficient. For a 
relative error of 1%, low-mass binaries are suitable up to red- 
shifts z ~ 0.1 — 1 while high-mass binaries can be observed 
up to z ~ 0.01 - 0.1. 
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ting bounds upon them. In figures 25p8 we present the 
maximal redshifts at which LISA can still measure the 
alternative theory parameters ^q, ^3/2 and ^ 2 f° r 
certain mass configurations with relative accuracies of 
a = 10% and a = 1%. Since for and vf/ 1/2 the fidu- 
cial 2PN phase coefficients are zero, we do not consider 
them. We checked that the error roughly scales with the 
redshift. For a relative accuracy of 1%, Vl^o is measurable 
up to redshifts of z ~ 1 — 10 for low-mass binaries and up 
to redshifts of z ~ 0.01 — 0.1 for high-mass binaries. 
^3/2 and ^2 can all be detected with a relative accuracy 
of 1% up to redshifts of z ~ 0.1 — 1 for low masses and 
z ~ 0.01 — 0.1 for high masses. For \&o> the use of the 
FWF improves the maximal redshifts by about a factor of 
2 for low masses and up to a factor of 10 for high masses, 
while the maximal redshifts are improved by almost an 
order of magnitude for the rest of the alternative theory 
parameters. If we were lucky and LISA could find a low- 
mass black hole binary at very low redshift z — 0.1, we 
would be able to recover the alternative theory parame- 
ters with ~ 10 times smaller errors than given in tables 

Ixvllxxl 



D. Example: Lower bound on graviton mass 

In order to compare our results with previous work in 
the field, we present here a lower bound on the graviton 
mass from our results at redshift z — 1. According to [5] , 
the effect of a massive graviton can be accounted for by 
introducing a gravitational wave phase correction 



A* 



MG 



= -pv 



-3/5^.-3/2 



(62) 



where x is the dimensionless frequency, v is the symmet- 
ric mass ratio and the parameter (3 is defined as 



A ^ 256^ /5A ^- 



(66) 



We take the fiducial value (3 = 0, thus the error A/3 
sets an upper bound on possible values for (3. A lower 
bound on the Compton wavelength of the graviton can 
then be calculated as 



\ g > 



'256 G ir 2 D{z)Mv 
~3~c2 (l + z)A*L 



(67) 



where M is the redshifted total mass of the binary. We 
find that optimal lower bounds on X g originate from a 
(3 x 10 6 + 1 x 10 7 )M Q binary for the FWF and from a 
(1 x 10 6 +1 x 10 6 ) M Q binary for the RWF. Including all six 
alternative theory parameters the resulting average 
bounds are X g > 1.2x 10 21 cm (FWF) and X g > 7.8x 10 20 
cm (RWF). These bounds are both lower than the one 
Yagi and Tanaka [TS] found (X g > 3.1 x 10 21 cm) using 
the RWF and simple precession at a distance of 3 Gpc; 
this is because the presence of the other five alternative 
theory parameters increases correlations among the pa- 
rameters. If we consider only one correction parameter 
which among other things accounts for massive grav- 
ity, the bounds increase to A g > 7.6 x 10 21 cm (FWF) 
and A s > 4.9 x 10 21 cm (RWF). The RWF bound is 
slightly higher than the one by Yagi and Tanaka for a 
(10 6 + 1O 7 )M0 binary; for this mass configuration we 
found a lower RWF bound of X g > 2.8 x 10 21 cm. Cor- 
nish et al. [2U] found a similar optimal RWF bound of 
A 9 > 3.8 x 10 21 cm. The use of the FWF improves the 
bound on the graviton Compton wavelength by a factor 
of ~ 1.5 with respect to the RWF, regardless whether 
only one or all the alternative theory parameters are in- 
cluded into the simulations. 



G ir 2 D{z)M 
J? A 2 (l + z) ' 



Here X g is the Compton wavelength of the graviton, z is 
the redshift, M = (1 + z)Mv 3/5 is the measured chirp 
mass affected by redshift, and D{z) is the distance given 
as 



£>(*) = (! + *)- 



dz' 



H Jo (1 + z') 2 \AM1 + ^') 3 + ^a' 

(64) 

where Ho, VLm and Vt\ are defined as in section [TTJ. In our 
framework, this is similar to the correction in Eq. (1401) : 



A* 



MG 



256^ 



Hence the errors on (3 and vf^ can be related with 



(63) VIII. CONCLUSION AND OUTLOOK 

We analyzed the expected measurement error distribu- 
tions of 17 different mass configurations of supermassive 
black hole binaries with masses between 10 5 — 10 8 M Q . 
We found that the black hole binaries can roughly be 
divided into two groups: low-mass binaries with M < 
10 7 M Q and high-mass binaries with M > 10 7 M Q . Com- 
paring the results of the simulations using the FWF and 
the RWF, we found that the RWF errors on the alter- 
native theory parameters ^ are a factor of ~ 100 times 
higher than the FWF errors for high-mass binaries, while 
they are almost comparable for low-mass binaries. Due 
to the dilution of the available information through the 
introduction of six extra parameters, the original param- 
eters lose accuracy. For masses and spins this is only a 

(65) 

factor of 1.2-5 for both low- and high-mass binaries re- 
gardless of whether the FWF or RWF is used. The loss 
of accuracy on the position of the black hole binary on 
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the sky is at maximum 10% for low-mass binaries and 
up to a factor of 5 for high-mass binaries. However, the 
accuracy of the luminosity distance is affected more se- 
riously for high-mass binaries, using the RWF results in 
a loss of a factor of ~ 50 — 1000 while using the FWF 
reduces it to factors of ~ 10 — 100. For low-mass binaries 
it is only about a factor of 2 worse. The use of the FWF 
is therefore mandatory for high-mass binaries, while the 
parameter estimation is more efficient for low-mass bina- 
ries and only up to a factor of 5 times worse when the 
RWF is used instead of the FWF. 

Since the error distributions were all calculated at fixed 
redshift z = 1 but the errors increase with redshift, we 
computed typical maximal redshifts up to which the al- 
ternative theory parameters are detectable with a relative 
accuracy smaller than 1% for the best 5% of the bina- 
ries in the sample. We found that for a deviation of 1% 
from the fiducial value, is detectable up to redshifts 
of z ~ 1 — 10 for low total masses and up to z ~ 0.01 — 0.1 
for high total masses. The rest of the alternative theory 
parameters ^3/2 and ^ with a fiducial 2PN phase 
coefficient unequal zero are detectable up to redshifts of 
z ~ 0.1 — 1 for low-mass binaries and z ~ 0.01 — 0.1 for 
high-mass binaries with the same relative accuracy. The 
use of the FWF improves the maximal redshifts up to a 
factor of 10 for high total masses. 

The FWF enables us to increase the optimal lower 
bound on the Compton wavelength of the graviton by 
about a factor of 1.5 compared to the one reached by 
the RWF. We achieve an optimal lower bound of X g > 
7.6 x 10 21 cm if only the alternative theory parameter ^1 
is considered. 

Since ESA is now considering a reduced LISA mission, 
certain aspects of the mission have to be reassessed. Sev- 
eral reduced variants of LISA are currently reviewed and 
will have an effect on the measurement accuracy of the 
mission [18]. It is thus of great importance to use as 
accurate waveform templates as possible to restore the 
lost accuracy with computational power on Earth. The 
use of the FWF improves the accuracy of the alterna- 
tive theory parameters by at least an order of magni- 
tude compared to the RWF. As shown by [15], the use 
of hybrid inspiral-merger-ringdown templates instead of 
inspiral-only templates improves the accuracy by an or- 
der of magnitude for the RWF; it would be interesting to 
find out how much such templates are improved when the 
FWF is used. The accuracy can further be enhanced by 
about an order of magnitude when considering combined 
observations instead of just extracting alternative theory 
parameters from individual black hole binaries [IB] . Also 
effects of eccentric orbits should be accounted for to make 
the model more realistic. 

Future work could include the introduction of ampli- 
tude corrections such as in [8J, since certain alternative 
theories are only visible in the gravitational wave am- 
plitude (e.g. Chern-Simons-modified gravity [37]). Also, 
the underlying mechanism of spin precession should be 
analyzed for effects originating from possible alternative 



theories. In this paper we neglected the energy loss of 
black hole binaries through unexpected physical effects 
such as further degrees of freedom in the propagation 
of gravitational waves arising from additional polariza- 
tions (e.g. longitudinal modes). It would be interesting 
to create a parameterized model for these effects. Also, 
since we studied a search for modifications at different 
PN orders at the same time, one could use the results of 
this work to investigate how the use of next-to-leading 
order modifications of GR could affect the determination 
of alternative theory parameters. 
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Appendix A: The 2.5PN and 3PN orbital frequency 
evolution equations 

The inclusion of dipole radiation corrections propor- 
tional to x^ 1 requires the knowledge of higher PN or- 
ders to be consistent to 2PN order, namely 2.5PN and 
3PN contributions. Since the current 2.5PN expansion 
just considers spin-orbit contributions and no spin-spin 
effects, and the 3PN expansion does not account for any 
spin coupling effects at all, these are of course only ap- 
proximations. 



1. 2.5PN from Blanchet et al. 2006 

Blanchet et al. 2006 [22 compute the angular fre- 
quency evolution for a binary with symmetric mass ratio 
v as 
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where u) — c 3 /(GM) x 3 ^ 2 , 5m = mi — m^ is the mass 
difference and m — mi + m 2 is the total mass. The spin 
interaction terms are expressed with 



18 



ft 

m 2 
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mi 



= 5 • J, £, = S ■ /, 

. (A2) 

where S = Si + S2 is the total spin and Z = j^y is the 
angular momentum unit vector. This enables us to write 
Eq. (Al) in the same form as Eq. (18), and we recover 
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2. 3PN without spin terms from Blanchet et al. 
2002 

In Luc Blanchet 's living review [21] (see also [38U40] ) . 
the 3PN expression for the total energy of non-spinning 
compact binaries can be found to be 
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Here v is the symmetric mass ratio and C — 0.577.. is 
the Euler constant. The logarithm in dE/dt will lead to 
a logarithmic term in the 3PN expansion. The PN coef- 
ficients bi can be recovered by computing the frequency 
evolution as a series in the dimensionless frequency x in 
the adiabatic approximation: 
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Appendix B: Tables 

I 



TABLE I: Median, 5% and 95% quantiles of the estimated measurement errors on m 1 for different mass configurations at 
rcdshift z = 1 with alternative theory parameters included. 
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TABLE II: 


Median, 5% and 95% quantiles of the estimated measurement 


errors on mi for different mass confij 


durations at 


redshift z = 


1 without considering alternative theory parameters. 
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TABLE III: Median, 5% and 95% quantiles of the estimated measurement errors on m 2 for different mass configurations at 
rcdshift 2 = 1 with alternative theory parameters included. 
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TABLE IV: Median, 5% and 95% quantiles of the estimated measurement errors on mi for different mass configurations at 
redshift z = 1 without considering alternative theory parameters. 
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3 x 10 5 


3 x 


10 4 


7.9 x 10" 5 


6.0 x 10~ b 


2.4 x 10" 4 


1.7 x 10" 4 


2.5 x 10" 3 


2.5 x 10" 3 


3 x 10 5 


1 X 


10 5 


1.5 x 10~ 4 


1.2 x 10~ 4 


6.7 x 10- 4 


3.9 x 10~ 4 


6.3 x 10" 3 


6.3 x 10" 3 


3 x 10 5 


3 x 


10 5 


1.8 x 10~ 4 


1.3 x 10~ 4 


1.2 x 10" 3 


2.8 x 10~ 4 


9.0 x 10" 3 


9.0 x 10~ 3 


1 x 10 6 


1 X 


10 5 


1.7 x 10~ 4 


1.1 x 10~ 4 


5.0 x 10" 4 


3.1 x 10" 4 


3.7 x 10" 3 


3.7 x 10" 3 


1 x 10 6 


3 x 


10 5 


3.0 x 10~ 4 


2.1 x 10~ 4 


1.1 x 10" 3 


6.5 x 10~ 4 


1.1 x 10~ 2 


1.1 x 10- 2 


1 x 10 6 


1 X 


10 6 


3.6 x 10~ 4 


2.6 x 10~ 4 


2.2 x 10" 3 


6.7 x 10~ 4 


2.1 x 10~ 2 


2.1 x 10~ 2 


3 x 10 6 


3 x 


10 5 


2.9 x 10~ 4 


1.6 x 10~ 4 


8.6 x 10" 4 


4.6 x 10" 4 


8.9 x 10" 3 


8.9 x 10" 3 


3 x 10 6 


1 X 


10 6 


7.0 x 10~ 4 


4.0 x 10~ 4 


2.9 x 10" 3 


9.2 x 10~ 4 


2.5 x 10~ 2 


2.5 x 10" 2 


3 x 10 6 


3 x 


10 6 


1.2 x 10" 3 


3.6 x 10- 4 


8.3 x 10" 3 


8.1 x 10" 4 


9.8 x 10- 2 


9.8 x 10" 2 


1 x 10 7 


1 X 


10 6 


1.1 x 10" 3 


3.9 x 10~ 4 


3.2 x 10" 3 


9.8 x 10~ 4 


3.3 x 10~ 2 


3.3 x 10- 2 


1 x 10 7 


3 x 


10 6 


2.8 x 10" 3 


7.2 x 10" 4 


1.3 x 10" 2 


1.5 x 10~ 3 


0.12 


0.12 


1 x 10 7 


1 X 


10 7 


2.3 x 10~ 2 


2.1 x 10" 3 


0.17 


5.6 x 10" 3 


1.7 


1.7 


3 x 10 7 


3 x 


10 6 


3.8 x 10~ 2 


2.6 x 10~ 3 


0.13 


5.2 x 10" 3 


0.56 


0.56 


3 x 10 7 


1 X 


10 7 


0.41 


5.8 x 10" 3 


1.4 


1.3 x 10" 2 


4.9 


4.9 


3 x 10 7 


3 x 


10 7 


3.1 


5.4 x 10~ 2 


21 


0.13 


250 


250 


1 x 10 s 


1 X 


10 7 


oo 


0.34 


oo 


1.0 


oo 


oo 


1 x 10 8 


3 x 


10 7 


oo 


3.0 


oo 


11 


oo 


oo 



TABLE V: Median, 5% and 95% quantiles of the estimated measurement errors on xi f° r different mass configurations at 
rcdshift z = 1 with alternative theory parameters included. 

m\[Mo\ m2[Mo\ Axi with corrections 

5%-quantile Median 95%-quantile 









RWF 




FWF 


RWF 




FWF 


RWF 


FWF 


3 x 10 5 


3 x 


10 4 


3.3 x 10" 


-4 


2.3 x 10~ 4 


1.1 x 10" 


-3 


7.1 x 10~ 4 


4.8 x 10" 3 


2.8 x 10" 3 


3 x 10 5 


1 X 


10 s 


8.8 x 10- 


-4 


6.0 x 10~ 4 


3.1 x 10- 


■3 


2.1 x 10~ 3 


1.5 x 10" 2 


1.0 x 10~ 2 


3 x 10 5 


3 x 


10 5 


1.9 x 10" 


-3 


1.5 x 10~ 3 


1.4 x 10- 


-2 


9.3 x 10" 3 


0.24 


0.12 


1 x 10 6 


1 X 


10 5 


6.8 x 10" 


-4 


4.3 x 10~ 4 


2.0 x 10- 


■3 


1.2 x 10" 3 


7.7 x 10" 3 


4.0 x 10~ 3 


1 x 10 6 


3 x 


10 s 


1.4 x 10- 


-3 


9.2 x 10~ 4 


4.7 x 10- 


■3 


3.0 x 10" 3 


2.6 x 10" 2 


1.2 x 10" 2 


1 x 10 6 


1 X 


10 6 


3.3 x 10- 


-3 


2.4 x 10" 3 


2.8 x 10- 


-2 


1.7 x 10~ 2 


0.58 


0.22 


3 x 10 6 


3 x 


10 5 


1.3 x 10- 


-3 


6.9 x 10- 4 


3.7 x 10- 


■3 


1.9 x 10" 3 


1.6 x 10" 2 


6.6 x 10" 3 


3 x 10 6 


1 X 


10 6 


3.0 x 10- 


-3 


1.4 x 10" 3 


1.3 x 10- 


-2 


4.9 x 10~ 3 


6.2 x 10~ 2 


2.4 x 10" 2 


3 x 10 6 


3 x 


10 6 


1.7 x 10- 


-2 


7.9 x 10" 3 


0.17 




5.4 x 10~ 2 


2.4 


0.58 


1 x 10 7 


1 X 


10 6 


5.3 x 10- 


-3 


1.5 x 10" 3 


1.8 x 10- 


-2 


4.3 x 10" 3 


7.4 x 10~ 2 


1.4 x 10~ 2 


1 x 10 7 


3 x 


10 6 


1.1 x 10- 


-2 


2.3 x 10" 3 


6.2 x 10- 


-2 


9.1 x 10~ 3 


0.33 


5.1 x 10" 2 


1 x 10 7 


1 X 


10 7 


0.62 




8.3 x 10~ 2 


4.1 




0.45 


26 


3.5 


3 x 10 7 


3 x 


10 6 


0.1 




6.6 x 10" 3 


0.41 




2.2 x 10~ 2 


2.1 


7.3 x 10~ 2 


3 x 10 7 


1 X 


10 7 


1.4 




2.1 x 10~ 2 


4.4 




9.0 x 10" 2 


16 


0.45 


3 x 10 7 


3 x 


10 7 


14 




3.1 


81 




13 


1.5 x 10 3 


85 


1 x 10 8 


1 X 


10 7 


oo 




0.76 


oo 




2.7 


oo 


10 


1 x 10 8 


3 x 


10 7 


oo 




8.8 


oo 




28 


oo 


270 



21 



TABLE VI: Median, 5% and 95% quantiles of the estimated measurement errors on xi for different mass configurations at 
redshift 2 = 1 without considering alternative theory parameters. 



mi[M ] 


m 2 [M Q ] 


RWF 


5%-quantile 

FWF 




Axi without corrections 

Median 
RWF FWF 




95%-quantile 
RWF FWF 


O X 1U 


3 x 


10 


1 q v in 
i.y x iu 


-4 


1.2 x 


10" 


-4 


4.3 x 10 


2.6 x 10" 


-4 


2.2 x 10 


2.2 X 10 


O X 1U 


1 X 


10 5 


A Q v in 


-4 


3.0 x 


10" 


-4 


1.5 x 10~ 3 


8.4 x 10" 


-4 


8.3 x 10 — 3 


8.3 x 10~ 3 


O X 1U 


3 x 


10 5 


1 n v in 


-3 


7.9 x 


10" 


-4 


6.2 x 10~ 3 


3.6 x 10- 


-3 


8.4 x 10 -2 


8.4 x 10~ 2 


1 v 1 f)6 

± X -LU 


1 X 


10 5 


q 7 v in 


-4 


1.9 x 


10" 


-4 


9.0 x 10~ 4 


4.6 x 10" 


-4 


2.6 x 10 -3 


2.6 x 10~ 3 


1 X 10 6 


3 x 


10 5 


8.5 x 10 


-4 


5.0 x 


io- 


-4 


2.2 x 10" 3 


1.2 x 10- 


-3 


1.4 x 10~ 2 


1.4 x 10~ 2 


1 X 10 6 


1 X 


10 6 


1.6 x 10 


-3 


1.2 x 


10" 


-3 


1.1 x 10~ 2 


6.0 x 10- 


-3 


0.17 


0.17 


3 x 10 6 


3 x 


10 s 


4.7 x 10 


-4 


2.2 x 


io- 


-4 


1.2 x 10" 3 


5.7 x 10- 


-4 


5.5 x 10~ 3 


5.5 x 10" 3 


3 x 10 6 


1 X 


10 6 


1.2 x 10 


-3 


6.2 x 


io- 


-4 


3.7 x 10" 3 


1.6 x 10- 


-3 


2.8 x 10~ 2 


2.8 x 10~ 2 


3 x 10 6 


3 x 


10 6 


4.2 x 10 


-3 


2.6 x 


io- 


-3 


4.1 x 10~ 2 


1.5 x 10- 


-2 


0.68 


0.68 


1 x 10 7 


1 X 


10 6 


1.4 x 10 


-3 


4.5 x 


io- 


-4 


3.4 x 10~ 3 


1.0 x 10- 


-3 


1.6 x 10~ 2 


1.6 x 10~ 2 


1 x 10 7 


3 x 


10 6 


3.6 x 10 


-3 


1.1 X 


io- 


-3 


1.5 x 10~ 2 


3.0 x 10- 


-3 


0.1 


0.1 


1 x 10 7 


1 X 


10 7 


0.13 




2.1 x 


io- 


-2 


1.1 


0.13 




9.7 


9.7 


3 x 10 7 


3 x 


10 6 


3.8 x 10 


-2 


2.7 x 


io- 


-3 


0.15 


6.8 x 10- 


-3 


0.69 


0.69 


3 x 10 7 


1 X 


10 7 


0.45 




9.2 x 


io- 


-3 


1.8 


3.2 x 10- 


-2 


6.7 


6.7 


3 x 10 7 


3 x 


10 7 


13 




1. 


1 




75 


5.1 




1.0 x 10 3 


1.0 x 10 3 


1 x 10 8 


1 X 


10 7 


oo 




0.25 




oo 


1.3 




oo 


oo 


1 x 10 8 


3 x 


10 7 


oo 




4.3 




oo 


16 




oo 


oo 



TABLE VII: Median, 5% and 95% quantiles of the estimated measurement errors on X2 for different mass configurations at 
redshift 2 = 1 with alternative theory parameters included. 



mi[M ] 


m 2 [M Q ] 






A%2 with corrections 










o/o- 


-quantile 




Median 


95%-quantile 








RWF 


FWF 


RWF 


FWF 


RWF 


FWF 


3 x 10 5 


3 x 


10 4 


2.0 x 10" 3 


1.2 X 10" 3 


1.5 x 10~ 2 


1.0 x 10" 2 


0.16 


8.9 x 10~ 2 


3 x 10 5 


1 X 


10 5 


1.7 x 10" 3 


1.1 x 10" 3 


1.0 x 10~ 2 


6.9 x 10~ 3 


8.4 x 10~ 2 


4.8 x 10~ 2 


3 x 10 5 


3 x 


10 s 


2.0 x 10" 3 


1.5 x 10~ 3 


1.4 x 10~ 2 


9.0 x 10~ 3 


0.23 


0.13 


1 x 10 6 


1 X 


10 5 


3.6 x 10" 3 


2.1 x 10~ 3 


2.6 x 10" 2 


1.6 x 10" 2 


0.26 


0.13 


1 x 10 6 


3 x 


10 s 


3.7 x 10" 3 


2.2 x 10~ 3 


1.9 x 10~ 2 


1.2 x 10" 2 


0.16 


7.3 x 10~ 2 


1 x 10 6 


1 X 


10 6 


3.6 x 10" 3 


2.5 x 10" 3 


2.8 x 10~ 2 


1.7 x 10~ 2 


0.54 


0.22 


3 x 10 6 


3 x 


10 5 


5.1 x 10" 3 


2.7 x 10" 3 


4.5 x 10" 2 


2.4 x 10" 2 


0.59 


0.21 


3 x 10 6 


1 X 


10 6 


6.0 x 10~ 3 


2.9 x 10" 3 


3.5 x 10~ 2 


1.7 x 10~ 2 


0.24 


0.1 


3 x 10 6 


3 x 


10 6 


1.6 x 10~ 2 


7.5 x 10" 3 


0.16 


5.6 x 10" 2 


2.6 


0.56 


1 x 10 7 


1 X 


10 6 


1.9 x 10~ 2 


4.8 x 10~ 3 


0.2 


5.7 x 10~ 2 


1.9 


0.37 


1 x 10 7 


3 x 


10 6 


2.2 x 10~ 2 


5.1 x 10~ 3 


0.15 


3.5 x 10~ 2 


1.3 


0.26 


1 x 10 7 


1 X 


10 7 


0.6 


7.8 x 10~ 2 


4.1 


0.47 


26 


3.3 


3 x 10 7 


3 x 


10 6 


0.39 


2.5 x 10~ 2 


3.7 


0.24 


37 


1.8 


3 x 10 7 


1 X 


10 7 


2.4 


3.7 x 10~ 2 


14 


0.3 


70 


1.7 


3 x 10 7 


3 x 


10 7 


17 


3.5 


84 


12 


1.3 x 10 3 


88 


1 x 10 s 


1 X 


10 7 


oo 


3.3 


oo 


25 


oo 


130 


1 x 10 8 


3 x 


10 7 


oo 


27 


oo 


130 


oo 


1.2 x 10 3 



TABLE VIII: Median, 5% and 95% quantiles of the estimated measurement errors on \2 for different mass configurations at 
redshift 2 = 1 without considering alternative theory parameters. 



mi[M ] 


m 2 [M Q ] 




A%2 without corrections 




5%-quantile 


Median 


95%-quantile 




^ „ A 


RWF FWF 

T^^L TT; 7\ Z — ~ — 7\ 


RWF FWF 

~ - TT^. 1 Z — - Z~ZZ 9 


RWF FWF 



3 x 10 5 


3 x 


10 4 


8.3 x 10- 4 


5.0 x 10" 4 


2.7 x IO" 3 


1.7 x 10~ 3 


1.9 x 10~ 2 


1.9 x 10" 2 


3 x 10 5 


1 X 


10 s 


7.4 x 10~ 4 


4.9 x 10~ 4 


3.4 x 10" 3 


2.0 x 10~ 3 


2.0 x 10~ 2 


2.0 x 10~ 2 


3 x 10 5 


3 x 


10 5 


9.0 x 10~ 4 


7.0 x 10~ 4 


6.6 x 10" 3 


3.7 x 10" 3 


7.6 x 10~ 2 


7.6 x 10~ 2 


1 x 10 6 


1 X 


10 5 


1.8 x 10" 3 


8.2 x 10~ 4 


5.6 x 10" 3 


3.0 x 10~ 3 


3.5 x 10" 2 


3.5 x 10~ 2 


1 x 10 6 


3 x 


10 s 


1.4 x 10" 3 


9.0 x 10~ 4 


4.8 x 10" 3 


2.8 x 10" 3 


3.5 x 10~ 2 


3.5 x 10~ 2 


1 x 10 6 


1 X 


10 6 


1.5 x 10" 3 


1.1 x 10" 3 


1.1 x 10" 2 


6.2 x 10" 3 


0.18 


0.18 


3 x 10 6 


3 x 


10 5 


2.2 x 10" 3 


1.1 x 10" 3 


7.4 x 10" 3 


3.8 x 10" 3 


5.9 x 10~ 2 


5.9 x 10~ 2 


3 x 10 6 


1 X 


10 6 


2.1 x 10" 3 


1.2 x 10" 3 


7.5 x 10" 3 


3.6 x 10" 3 


4.8 x 10" 2 


4.8 x 10~ 2 


3 x 10 6 


3 x 


10 6 


4.6 x 10" 3 


2.5 x 10~ 3 


4.1 x 10~ 2 


1.5 x 10" 2 


0.68 


0.68 


1 x 10 7 


1 X 


10 6 


6.4 x 10" 3 


2.0 x 10~ 3 


3.4 x 10~ 2 


9.4 x 10" 3 


0.19 


0.19 


1 x 10 7 


3 x 


10 6 


7.4 x 10" 3 


2.4 x 10" 3 


3.2 x 10" 2 


9.3 x 10~ 3 


0.18 


0.18 


1 x 10 7 


1 X 


10 7 


0.12 


2.2 x 10" 2 


1.1 


0.13 


9.8 


9.8 


3 x 10 7 


3 x 


10 6 


0.16 


1.0 x 10" 2 


1.4 


6.5 x 10~ 2 


7.0 


7.0 


3 x 10 7 


1 X 


10 7 


0.53 


1.7 x 10~ 2 


4.3 


9.8 x 10" 2 


23 


23 


3 x 10 7 


3 x 


10 7 


16 


1.4 


70 


4.3 


910 


910 


1 x 10 8 


1 X 


10 7 


oo 


1.5 


oo 


12 


oo 


oo 


1 x 10 8 


3 x 


10 7 


oo 


11 


oo 


53 


oo 


oo 



22 



TABLE IX: Median, 5% and 95% quantiles of the estimated measurement errors on 2a for different mass configurations at 
rcdshift z = 1 with alternative theory parameters included. 



mi[M ] 


m 2 [M ] 


5%-quantile 
RWF FWF 


2ct[ ; ] with corrections 
Median 
RWF FWF 


95%-quantile 
RWF FWF 


3 x 10 & 


3 x 10 4 


7.8 


4.8 


24 


16 


92 


74 


3 x 10 5 


1 x 10 5 


5.1 


3.2 


23 


15 


100 


86 


3 x 10 5 


3 x 10 5 


6.9 


4.5 


32 


24 


130 


110 


1 x 10 6 


1 x 10 5 


\\ 


7.9 


35 


20 


140 


92 


1 x 10 6 


3 x 10 5 


11 


7^3 


36 


24 


140 


100 


1 x 10 6 


1 x 10 6 


11 


7.3 


46 


31 


190 


150 


3 x 10 6 


3 x 10 5 


13 


7.3 


41 


22 


140 


87 


3 x 10 6 


1 x 10 6 


13 


6.6 


47 


25 


220 


130 


3 x 10 6 


3 x 10 6 


15 


7.6 


77 


41 


400 


229 


1 x 10 7 


1 x 10 6 


28 


9.1 


81 


26 


459 


130 


1 x 10 7 


3 x 10 6 


24 


8.4 


120 


41 


800 


200 


1 x 10 7 


1 x 10 7 


57 


18 


440 


120 


3.6 x 10 3 


750 


3 x 10 7 


3 x 10 6 


170 


25 


630 


75 


6.8 x 10 3 


310 


3 x 10 7 


1 x 10 7 


409 


38 


4.4 x 10 3 


190 


6.3 x 10 4 


1.2 x 10 3 


3 x 10 7 


3 x 10 7 


1.1 x 10 4 


590 


4.0 x 10 5 


4.9 x 10 3 


1.3 x 10 7 


3.0 x 10 4 


1 x 10 8 


1 x 10 7 


oo 


2.2 x 10 3 


oo 


8.2 x 10 3 


oo 


4.3 x 10 4 


1 x 10 s 


3 x 10 7 


oo 


1.6 x 10 4 


oo 


8.3 x 10 4 


oo 


7.0 x 10 5 



TABLE X: Median, 5% and 95% quantiles of the estimated measurement errors on 2a for different mass configurations at 
redshift z = 1 without considering alternative theory parameters. 



mi [M ] 


m 2 [M e ] 


5%-quantile 
RWF FWF 


2a['] without corrections 

Median 
RWF FWF 


95%-quantile 
RWF FWF 


3 x 10 b 


3 x 10 4 


7.4 


4.6 


22 


14 


86 


86 


3 x 10 5 


1 x 10 5 


4.7 


2.9 


20 


13 


97 


97 


3 x 10 5 


3 x 10 5 


6.0 


3.9 


27 


20 


110 


110 


1 x 10 6 


1 x 10 5 


7.5 


6.7 


32 


18 


130 


130 


1 x 10 6 


3 x 10 5 


9.0 


6.7 


33 


22 


130 


130 


1 x 10 6 


1 x 10 6 


8.6 


6.2 


40 


26 


160 


160 


3 x 10 6 


3 x 10 5 


7.2 


5.0 


36 


20 


130 


130 


3 x 10 6 


1 x 10 6 


11 


5.8 


38 


21 


160 


160 


3 x 10 6 


3 x 10 6 


10 


5.3 


53 


29 


280 


280 


1 x 10 7 


1 x 10 6 


18 


8.2 


61 


23 


280 


280 


1 x 10 7 


3 x 10 6 


20 


7.5 


89 


31 


420 


420 


1 x 10 7 


1 x 10 7 


33 


13 


200 


76 


1.4 x 10 3 


1.4 x 10 3 


3 x 10 7 


3 x 10 6 


100 


21 


370 


60 


2.8 x 10 3 


2.8 x 10 3 


3 x 10 7 


1 x 10 7 


160 


27 


1.3 x 10 3 


130 


1.3 x 10 4 


1.3 x 10 4 


3 x 10 7 


3 x 10 7 


3.9 x 10 3 


420 


1.6 x 10 5 


2.7 x 10 3 


2.3 x 10 6 


2.3 x 10 6 


1 x 10 s 


1 x 10 7 


oo 


1.2 X 10 3 


oo 


4.5 x 10 3 


oo 


oo 


1 x 10 8 


3 x 10 7 


oo 


7.2 x 10 3 


oo 


3.5 x 10 4 


oo 


oo 



TABLE XI: Median, 5% and 95% quantiles of the estimated measurement errors on 26 for different mass configurations at 
rcdshift 2 = 1 with alternative theory parameters included. 



mi [M Q ] 


m 2 [M ] 






26['] 


with corrections 






5%-quantile 


Median 


95%-quantile 






RWF 


FWF 


RWF 


FWF 


RWF 


FWF 


3 x 10 5 


3 x 10 4 


1.0 


0.59 


4.8 


2.8 


14 


7.5 


3 x 10 5 


1 x 10 5 


0.86 


0.49 


4.2 


2.4 


20 


11 


3 x 10 5 


3 x 10 5 


1.3 


0.77 


6.3 


3.9 


24 


17 


1 x 10 6 


1 x 10 5 


2.3 


1.0 


11 


5.0 


26 


11 


1 x 10 6 


3 x 10 5 


2.1 


1.2 


9.1 


4.9 


28 


17 


1 x 10 6 


1 x 10 6 


1.9 


1.1 


9.2 


5.6 


33 


23 


3 x 10 6 


3 x 10 5 


2.0 


0.93 


10 


4.9 


30 


14 


3 x 10 6 


1 x 10 6 


2.0 


0.95 


9.3 


4.5 


33 


16 


3 x 10 6 


3 x 10 6 


2.2 


1.0 


12 


5.6 


52 


27 


1 x 10 7 


1 x 10 6 


3.1 


1.1 


17 


5.4 


72 


18 


1 x 10 7 


3 x 10 6 


4.0 


1.4 


18 


6.1 


83 


25 


1 x 10 7 


1 x 10 7 


13 


4.0 


60 


17 


320 


89 


3 x 10 7 


3 x 10 6 


33 


4.2 


120 


15 


810 


55 


3 x 10 7 


1 x 10 7 


82 


9.3 


520 


29 


5.7 x 10 3 


130 


3 x 10 7 


3 x 10 7 


1.1 x 10 3 


98 


3.0 x 10 


4 560 


8.9 x 10 5 


2.6 x 10 3 


1 x 10 8 


1 x 10 7 


oo 


459 


oo 


1.7 x 10 3 


oo 


7.1 x 10 3 


1 x 10 s 


3 x 10 7 


oo 


2.8 x 10 3 


oo 


1.3 x 10 4 


oo 


1.0 x 10 5 



23 



TABLE XII 


: Median, 5% and 95% quantiles of the estimated measurement 


errors on 2b for different mass 


configurations at 


redshift z = 


1 without considering alternative theory parameters. 




















2M'l without corrections 








5%-quantile 




Median 


95%-quantilc 








RWF 


FWF 


KWr 




FWF 


RWF 


FWF 


3 x 10 5 




3 x 10 4 


u.y i 


u.oo 


a a 
4.0 




9 G. 
Z.O 


1 ^ 
10 




3 x 10 5 




1 x 10 5 


A 1 A 
U. i 4 


n A q 


3.9 




O 9 

z.z 


1 V 

1 ( 


1 V 

1 ( 


3 x 10 5 




3 x 10 5 


1 . 1 


u.oo 


5.6 




Q 9 
O.Z 


9fl 

zu 


9fl 


1 x 10 6 




1 x 10 5 


9 1 
z. 1 


U.y4 


10 




4.0 


9£ 

zu 


9 K 
ZO 


1 x 10 6 




3 x 10 5 


2.0 


1.0 


8.7 




4.6 


25 


25 


1 x 10 6 




1 x 10 6 


1.7 


0.96 


8.4 




4.9 


29 


29 


3 x 10 6 




3 x 10 5 


1.8 


0.85 


9.8 




4.6 


26 


26 


3 x 10 6 




1 x 10 6 


1.6 


0.85 


8.4 




4.0 


28 


28 


3 x 10 6 




3 x 10 6 


1.5 


0.77 


9.5 




4.6 


39 


39 


1 x 10 7 




1 x 10 6 


2.7 


0.87 


14 




4.8 


45 


45 


1 x 10 7 




3 x 10 6 


3.1 


1.1 


14 




4.9 


57 


57 


1 x 10 7 




1 x 10 7 


6.2 


2.1 


26 




9.8 


120 


120 


3 x 10 7 




3 x 10 6 


16 


2.7 


65 




12 


290 


290 


3 x 10 7 




1 x 10 7 


40 


4.2 


150 




16 


910 


910 


3 x 10 7 




3 x 10 7 


540 


49 


4.8 x 10 3 




260 


1.6 x 10 5 


1.6 x 10 5 


1 x 10 s 




1 x 10 7 


oo 


260 


oo 




900 


oo 


oo 


1 x 10 s 




3 x 10 7 


oo 


1.4 x 10 3 


oo 




4.9 x 10 3 


oo 


oo 


TABLE XIII: Median, 5% and 95% quantiles of the estimated measurement 


errors on <Il for different mass 


configurations at 


redshift z = 


1 with alternative theory parameters included. 












mi [M Q ] 


m 2 [M ] 






Ad L /d L 


with corrections 










5%-quantile 




Median 


95%-quantile 








RWF 


FWF 


RWF 




FWF 


RWF 


FWF 


3 x 10 b 


3 


x 10 4 


2.6 X 10--* 


1.6 x 10~ 3 


5.0 x 10" 3 




3.1 x 10~ 3 


1.5 X 10" 2 


9.2 x 10" 3 


3 x 10 5 


1 


x 10 5 


2.2 x 10~ 3 


1.5 x 10~ 3 


5.1 x 10~ 3 




3.4 x 10~ 3 


1.7 x 10~ 2 


1.2 x 10~ 2 


3 x 10 5 


3 


x 10 5 


2.6 x 10~ 3 


1.8 x 10~ 3 


7.4 x 10~ 3 




5.3 x 10~ 3 


2.2 x 10~ 2 


1.8 x 10~ 2 


1 x 10 6 


1 


x 10 5 


3.8 x 10~ 3 


2.4 x 10~ 3 


7.9 x 10~ 3 




5.1 x 10~ 3 


2.3 x 10~ 2 


1.3 x 10~ 2 


1 x 10 6 


3 


x 10 5 


3.2 x 10~ 3 


2.5 x 10- 3 


7.3 x 10- 3 




5.0 x 10~ 3 


2.5 x 10- 2 


1.6 x 10- 2 


1 x 10 6 


1 


x 10 6 


3.8 x 10~ 3 


2.5 x 10~ 3 


1.0 x 10~ 2 




7.3 x 10~ 3 


3.0 x 10~ 2 


2.3 x 10~ 2 


3 x 10 6 


3 


x 10 5 


4.6 x 10~ 3 


3.0 x 10~ 3 


9.8 x 10~ 3 




6.1 x 10~ 3 


2.9 x 10~ 2 


1.7 x 10~ 2 


3 x 10 6 


1 


x 10 6 


5.2 x 10~ 3 


3.4 x 10~ 3 


1.2 x 10~ 2 




7.4 x 10~ 3 


3.5 x 10~ 2 


2.0 x 10" 2 


3 x 10 6 


3 


x 10 6 


7.5 x 10~ 3 


4.6 x 10~ 3 


2.3 x 10~ 2 




1.3 x 10~ 2 


8.0 x 10~ 2 


4.2 x 10~ 2 


1 x 10 7 


1 


x 10 6 


1.4 x 10~ 2 


5.9 x 10~ 3 


3.6 x 10~ 2 




1.4 x 10~ 2 


0.12 


3.2 x 10" 2 


1 x 10 7 


3 


x 10 6 


2.6 x 10~ 2 


7.9 x 10~ 3 


6.9 x 10~ 2 




1.7 x 10~ 2 


0.22 


4.5 x 10~ 2 


1 x 10 7 


1 


x 10 7 


0.23 


4.4 x 10~ 2 


0.7 




9.3 x 10~ 2 


3.2 


0.25 


3 x 10 7 


3 


x 10 6 


0.97 


3.8 x 10~ 2 


4.5 




8.0 x 10~ 2 


23 


0.19 


3 x 10 7 


1 


x 10 7 


8.5 


0.11 


40 




0.27 


200 


0.71 


3 x 10 7 


3 


x 10 7 


630 


5.1 


1.7 x 10 4 




12 


2.9 x 10 5 


52 


1 x 10 s 


1 


x 10 7 


oo 


17 


oo 




73 


oo 


440 


1 x 10 8 


3 


x 10 7 


oo 


170 


oo 




819 


oo 


1.7 x 10 4 


TABLE XIV: Median, 5% and 95% quantiles of the estimated measurement 


errors on for different mass 


configurations at 


redshift z = 


1 without considering alternative theory paramefc 


ers. 










mi[M Q ] 


m 2 \M Q ] 






Adj^/d]^ without corrections 










5%-quantile 




Median 


95%-quantile 








RWF 


FWF 


RWF 




FWF 


RWF 


FWF 


3 x 10 5 


3 


x 10 4 


1.2 x 10" 3 


8.1 x 10~ 4 


2.7 x 10~ 3 




1.6 x 10" 3 


1.2 x 10~ 2 


1.2 x 10~ 2 


3 x 10 5 


1 


x 10 5 


8.7 x 10- 4 


5.6 x 10~ 4 


2.8 x 10~ 3 




1.7 x 10~ 3 


1.4 x 10~ 2 


1.4 x 10~ 2 


3 x 10 5 


3 


x 10 5 


1.2 x 10~ 3 


7.3 x 10~ 4 


4.5 x 10~ 3 




2.9 x 10~ 3 


1.7 x 10~ 2 


1.7 x 10~ 2 


1 x 10 6 


1 


x 10 5 


1.9 x 10~ 3 


9.5 x 10~ 4 


4.9 x 10~ 3 




2.4 x 10~ 3 


1.9 x 10~ 2 


1.9 x 10~ 2 


1 x 10 6 


3 


x 10 5 


1.8 x 10~ 3 


1.1 x 10~ 3 


5.0 x 10~ 3 




2.9 x 10~ 3 


2.1 x 10~ 2 


2.1 x 10" 2 


1 x 10 6 


1 


x 10 6 


1.8 x 10~ 3 


1.2 x 10~ 3 


6.2 x 10~ 3 




4.0 x 10~ 3 


2.4 x 10~ 2 


2.4 x 10~ 2 


3 x 10 6 


3 


x 10 5 


1.8 x 10~ 3 


1.0 x 10~ 3 


5.1 x 10~ 3 




2.6 x 10~ 3 


2.1 x 10" 2 


2.1 x 10" 2 


3 x 10 6 


1 


x 10 6 


2.0 x 10~ 3 


1.1 x 10~ 3 


5.8 x 10~ 3 




3.2 x 10~ 3 


2.2 x 10~ 2 


2.2 x 10~ 2 


3 x 10 6 


3 


x 10 6 


2.1 x 10~ 3 


1.4 x 10~ 3 


8.3 x 10~ 3 




4.7 x 10~ 3 


4.6 x 10" 2 


4.6 x 10~ 2 


1 x 10 7 


1 


x 10 6 


2.8 x 10~ 3 


1.3 x 10~ 3 


8.1 x 10~ 3 




3.3 x 10~ 3 


3.9 x 10~ 2 


3.9 x 10" 2 


1 x 10 7 


3 


x 10 6 


4.1 x 10~ 3 


2.2 x 10~ 3 


1.1 x 10~ 2 




4.7 x 10~ 3 


6.6 x 10~ 2 


6.6 x 10~ 2 


1 x 10 7 


1 


x 10 7 


1.4 x 10~ 2 


7.3 x 10~ 3 


6.1 x 10~ 2 




2.1 x 10~ 2 


0.27 


0.27 


3 x 10 7 


3 


x 10 6 


3.0 x 10~ 2 


5.0 x 10~ 3 


7.9 x 10~ 2 




1.1 x 10~ 2 


0.48 


0.48 


3 x 10 7 


1 


x 10 7 


0.21 


1.6 x 10~ 2 


0.75 




3.3 x 10~ 2 


3.1 


3.1 


3 x 10 7 


3 


x 10 7 


3.0 


0.23 


30 




0.56 


390 


390 


1 x 10 8 


1 


x 10 7 


oo 


0.41 


oo 




1.1 


oo 


oo 


1 x 10 8 


3 


x 10 7 


oo 


3.4 


oo 




9.9 


oo 


oo 



24 



TABLE XV: Median, 5% and 95% quantiles of the estimated measurement errors on for different mass configurations at 
redshift 2 = 1. 

mi[M ] m 2 [M Q ] A*_i 

5%-quantile Median 95%-quantile 











RWF 




FWF 


RWF 


FWF 


RWF 


FWF 


3 x 10 b 


3 


X 


10 4 


3.0 x 10" 


-6 


1.6 x 10~ b 


7.2 x 10~ b 


5.3 x 10~ b 


6.5 x 10 _t> 


5.0 x 10 _t> 


3 x 10 5 


1 


X 


10 5 


2.9 x 10" 


■6 


1.8 x 10~ 6 


7.7 x 10~ 6 


5.8 x 10~ 6 


6.1 x 10~ 5 


5.2 x 10~ 5 


3 x 10 5 


3 


X 


10 5 


5.3 x 10" 


■6 


3.1 x 10~ 6 


1.6 x 10~ 5 


1.1 x 10~ 5 


1.8 x 10~ 4 


1.3 x 10~ 4 


1 x 10 6 


1 


X 


10 5 


1.5 x 10" 


•5 


6.9 x 10 -6 


3.5 x 10~ 5 


2.2 x 10~ 5 


3.3 x 10~ 4 


2.4 x 10~ 4 


1 x 10 6 


3 


X 


10 5 


1.5 x 10" 


■5 


7.8 x 10~ 6 


3.5 x 10" B 


2.6 x 10~ 5 


2.7 x 10~ 4 


1.8 x 10~ 4 


1 x 10 6 


1 


X 


10 6 


5.1 x 10" 


-5 


2.2 x 10" B 


1.4 x 10~ 4 


8.7 x 10~ 5 


7.7 x 10~ 4 


5.0 x 10~ 4 


3 x 10 6 


3 


X 


10 s 


9.8 x 10" 


■5 


2.5 x 10~ 5 


2.0 x 10~ 4 


1.1 x 10~ 4 


1.2 x 10" 3 


6.4 x 10~ 4 


3 x 10 6 


1 


X 


10 6 


2.0 x 10" 


-4 


6.6 x 10~ 5 


4.0 x 10~ 4 


1.7 x 10~ 4 


1.1 x 10~ 3 


6.4 x 10~ 4 


3 x 10 6 


3 


X 


10 6 


8.8 x 10" 


-4 


2.0 x 10~ 4 


2.6 x 10~ 3 


7.5 x 10~ 4 


8.7 x 10~ 3 


2.4 x 10" 3 


1 x 10 7 


1 


X 


10 6 


2.6 x 10" 


-3 


3.1 x 10~ 4 


6.3 x 10" 3 


9.0 x 10~ 4 


1.9 x 10~ 2 


3.0 x 10" 3 


1 x 10 7 


3 


X 


10 6 


6.9 x 10" 


-3 


1.1 x 10" 3 


1.9 x 10~ 2 


2.0 x 10" 3 


5.6 x 10~ 2 


4.8 x 10" 3 


1 x 10 7 


1 


X 


10 7 


6.2 x 10" 


-2 


9.8 x 10" 3 


0.18 


2.5 x 10~ 2 


0.77 


6.0 x 10~ 2 


3 x 10 7 


3 


X 


10 6 


0.28 




8.9 x 10" 3 


1.4 


1.7 x 10~ 2 


8.1 


5.1 x 10~ 2 


3 x 10 7 


1 


X 


10 7 


2.3 




3.6 x 10~ 2 


14 


0.1 


71 


0.27 


3 x 10 7 


3 


X 


10 7 


190 




1.7 


5.9 x 10 3 


4.1 


1.1 x 10 B 


17 


1 x 10 8 


1 


X 


10 7 


oo 




4.6 


oo 


22 


oo 


170 


1 x 10 8 


3 


X 


10 7 


oo 




55 


oo 


280 


oo 


6.1 x 10 3 



TABLE XVI: Median, 5% and 95% quantiles of the estimated measurement errors on for different mass configurations at 
redshift 2 = 1. 

mi[M©] m2[M Q ] A^o 

5%-quantile Median 95%-quantile 











RWF 


FWF 


RWF 


FWF 


RWF 


FWF 


3 x 10 5 


3 


X 


10 4 


1.5 x 10" 3 


7.4 x 10~ 4 


3.2 x 10" 3 


2.1 x 10~ 3 


1.7 x 10~ 2 


1.2 x 10" 2 


3 x 10 5 


1 


X 


10 5 


1.3 x 10" 3 


7.2 x 10~ 4 


3.3 x 10" 3 


2.3 x 10" 3 


1.7 x 10~ 2 


1.4 x 10~ 2 


3 x 10 5 


3 


X 


10 5 


2.3 x 10" 3 


1.4 x 10" 3 


6.1 x 10~ 3 


4.5 x 10" 3 


4.3 x 10~ 2 


3.2 x 10~ 2 


1 x 10 6 


1 


X 


10 5 


4.8 x 10" 3 


1.9 x 10" 3 


1.1 x 10" 2 


6.3 x 10" 3 


7.0 x 10~ 2 


4.7 x 10~ 2 


1 x 10 6 


3 


X 


10 s 


4.9 x 10" 3 


2.4 x 10" 3 


1.1 x 10~ 2 


7.8 x 10" 3 


6.1 x 10" 2 


4.1 x 10~ 2 


1 x 10 6 


1 


X 


10 6 


1.4 x 10~ 2 


5.2 x 10~ 3 


3.8 x 10~ 2 


2.3 x 10" 2 


0.16 


0.11 


3 x 10 6 


3 


X 


10 5 


2.1 x 10~ 2 


4.4 x 10~ 3 


4.5 x 10" 2 


2.4 x 10" 2 


0.21 


0.11 


3 x 10 6 


1 


X 


10 6 


3.8 x 10~ 2 


1.1 x 10~ 2 


8.4 x 10~ 2 


3.6 x 10~ 2 


0.21 


0.12 


3 x 10 6 


3 


X 


10 6 


0.13 


2.2 x 10~ 2 


0.38 


0.12 


1.3 


0.37 


1 x 10 7 


1 


X 


10 6 


0.3 


3.5 x 10~ 2 


0.75 


0.14 


2.4 


0.43 


1 x 10 7 


3 


X 


10 6 


0.7 


0.12 


2.0 


0.27 


6.0 


0.67 


1 x 10 7 


1 


X 


10 7 


4.9 


0.96 


14 


2.4 


60 


6.2 


3 x 10 7 


3 


X 


10 6 


9.8 


0.78 


76 


1.8 


590 


4.9 


3 x 10 7 


1 


X 


10 7 


88 


2.7 


720 


8.9 


4.3 x 10 3 


25 


3 x 10 7 


3 


X 


10 7 


3.8 x 10 3 


120 


1.7 x 10 B 


300 


6.2 x 10 6 


1.5 x 10 3 


1 x 10 s 


1 


X 


10 7 


oo 


150 


oo 


980 


oo 


1.1 x 10 4 


1 x 10 8 


3 


X 


10 7 


oo 


2.1 x 10 3 


oo 


1.4 x 10 4 


oo 


4.1 x 10 5 



TABLE XVII: Median, 5% and 95% quantiles of the estimated measurement errors on *i/ 2 for different mass configurations 
at redshift 2 = 1. 

mi[M ] m 2 [M Q ] At 1/2 

5%-quantile Median 95%-quantile 







RWF 


FWF 


RWF 


FWF 


RWF 


FWF 


3 x 10 5 


3 x 10 4 


1.6 x 10" 2 


5.8 x 10" 3 


3.7 x 10" 2 


2.0 x 10" 2 


0.14 


0.11 


3 x 10 5 


1 x 10 5 


1.3 x 10~ 2 


7.0 x 10" 3 


3.4 x 10~ 2 


2.3 x 10~ 2 


0.14 


0.11 


3 x 10 5 


3 x 10 5 


1.8 x lO" 2 


1.1 x 10~ 2 


5.3 x 10~ 2 


4.0 x 10~ 2 


0.29 


0.21 


1 x 10 6 


1 x 10 5 


4.0 x 10~ 2 


1.5 x 10~ 2 


9.1 x 10" 2 


5.0 x 10~ 2 


0.45 


0.3 


1 x 10 6 


3 x 10 5 


3.7 x 10" 2 


1.9 x 10~ 2 


8.7 x 10~ 2 


6.2 x 10" 2 


0.42 


0.28 


1 x 10 6 


1 x 10 6 


9.8 x 10~ 2 


3.7 x 10" 2 


0.27 


0.16 


0.98 


0.7 


3 x 10 6 


3 x 10 5 


0.12 


3.1 x 10~ 2 


0.29 


0.15 


1.2 


0.65 


3 x 10 6 


1 x 10 6 


0.23 


7.2 x 10~ 2 


0.51 


0.23 


1.3 


0.73 


3 x 10 6 


3 x 10 6 


0.66 


0.12 


2.0 


0.64 


6.6 


2.1 


1 x 10 7 


1 x 10 6 


1.3 


0.17 


3.6 


0.75 


12 


2.4 


1 x 10 7 


3 x 10 6 


3.0 


0.55 


9.0 


1.3 


29 


3.5 


1 x 10 7 


1 x 10 7 


19 


4.1 


55 


10 


240 


29 


3 x 10 7 


3 x 10 6 


23 


3.1 


250 


7.8 


2.1 X 10 3 


23 


3 x 10 7 


1 x 10 7 


150 


10 


2.3 x 10 3 


37 


1.4 x 10 4 


110 


3 x 10 7 


3 x 10 7 


1.2 x 10 4 


430 


3.4 x 10 5 


1.1 x 10 3 


7.7 x 10 6 


6.4 x 10 3 


1 x 10 s 


1 x 10 7 


oo 


320 


oo 


2.9 x 10 3 


oo 


4.1 x 10 4 


1 x 10 8 


3 x 10 7 


oo 


4.5 x 10 3 


oo 


4.5 x 10 4 


oo 


1.4 x 10 6 
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TABLE XVIII: Median, 5% and 95% quantiles of the estimated measurement errors on for different mass configurations 
at redshift 2 = 1. 

mi[M Q ] m 2 [M Q ] A*i 

5%-quantile Median 95%-quantiie 







RWF 


FWF 


RWF 


FWF 


RWF 


FWF 


3 x 10 b 


3 x 10 4 


0.15 


3.5 x 10" 2 


0.33 


0.17 


0.92 


0.64 


3 x 10 5 


1 x 10 5 


0.14 


4.7 x 10 -2 


0.34 


0.2 


0.95 


0.79 


3 x 10 5 


3 x 10 5 


0.19 


6.8 x 10~ 2 


0.49 


0.29 


1.7 


1.3 


1 x 10 6 


1 x 10 5 


0.3 


7.9 x 10 -2 


0.65 


0.34 


2.1 


1.3 


1 x 10 6 


3 x 10 5 


0.29 


0.11 


0.64 


0.41 


1.9 


1.4 


1 x 10 6 


1 x 10 6 


0.47 


0.13 


1.2 


0.65 


3.7 


2.8 


3 x 10 6 


3 x 10 5 


0.61 


0.13 


1.4 


0.7 


4.1 


2.2 


3 x 10 6 


1 x 10 6 


0.78 


0.21 


1.8 


0.86 


4.7 


2.5 


3 x 10 6 


3 x 10 6 


1.3 


0.29 


4.0 


1.4 


16 


6.1 


1 x 10 7 


1 x 10 6 


3.0 


0.36 


8.8 


2.3 


30 


5.8 


1 x 10 7 


3 x 10 6 


5.1 


0.92 


14 


2.6 


45 


7.2 


1 x 10 7 


1 x 10 7 


26 


5.7 


76 


13 


330 


35 


3 x 10 7 


3 x 10 6 


58 


6.6 


570 


16 


4.0 x 10 3 


38 


3 x 10 7 


1 x 10 7 


409 


15 


3.3 x 10 3 


46 


1.9 x 10 4 


120 


3 x 10 7 


3 x 10 7 


2.1 x 10 4 


459 


4.3 x 10 5 


1.2 x 10 3 


2.1 x 10 7 


5.0 x 10 3 


1 x 10 s 


1 x 10 7 


oo 


919 


oo 


6.7 x 10 3 


oo 


6.8 x 10 4 


1 x 10 8 


3 x 10 7 


oo 


7.8 x 10 3 


oo 


5.9 x 10 4 


oo 


1.4 x 10 6 



TABLE XIX: Median, 5% and 95% quantiles of the estimated measurement errors on ^ 3 / 2 for different mass configurations 
at redshift 2 = 1. 

mi[M s ] m 2 [M©] A<£ 3/2 

5%-quantile Median 95%-quantile 







RWF 


FWF 


RWF 


FWF 


RWF 


FWF 


3 x 10 5 


3 x 10 4 


1.1 


0.16 


2.5 


1.5 


8.0 


5.7 


3 x 10 5 


1 x 10 5 


1.3 


0.18 


3.1 


1.8 


9.1 


7.3 


3 x 10 5 


3 x 10 5 


2.0 


0.3 


5.3 


2.9 


18 


14 


1 x 10 6 


1 x 10 5 


2.9 


0.33 


6.3 


3.5 


25 


17 


1 x 10 6 


3 x 10 5 


3.2 


0.44 


7.1 


4.5 


24 


16 


1 x 10 6 


1 x 10 6 


6.7 


0.85 


17 


8.6 


44 


32 


3 x 10 6 


3 x 10 5 


7.6 


0.48 


17 


9.1 


54 


28 


3 x 10 6 


1 x 10 6 


11 


1.4 


26 


12 


67 


35 


3 x 10 6 


3 x 10 6 


20 


1.7 


54 


19 


170 


71 


1 x 10 7 


1 x 10 6 


32 


2.7 


110 


32 


360 


81 


1 x 10 7 


3 x 10 6 


56 


11 


180 


37 


590 


100 


1 x 10 7 


1 x 10 7 


250 


61 


770 


150 


3.6 x 10 3 


459 


3 x 10 7 


3 x 10 6 


459 


71 


4.0 x 10 3 


190 


3.4 x 10 4 


490 


3 x 10 7 


1 x 10 7 


2.9 x 10 3 


140 


2.4 x 10 4 


490 


1.4 x 10 5 


1.5 x 10 3 


3 x 10 7 


3 x 10 7 


1.3 x 10 5 


3.2 x 10 3 


2.4 x 10 6 


1.1 x 10 4 


6.4 x 10 7 


6.8 x 10 4 


1 x 10 8 


1 x 10 7 


oo 


8.4 x 10 3 


oo 


4.4 x 10 4 


oo 


5.9 x 10 5 


1 x 10 s 


3 x 10 7 


oo 


6.0 x 10 4 


oo 


4.5 x 10 5 


oo 


1.4 x 10 7 



TABLE XX: Median, 5% and 95% quantiles of the estimated measurement errors on ^2 for different mass configurations at 
redshift 2 = 1. 

mi[M ] m 2 [M e ] A* 2 

5%-quantile Median 95%-quantile 







RWF 


FWF 


RWF 


FWF 


RWF 


FWF 


3 x 10 5 


3 x 10 4 


4.0 


0.36 


9.1 


5.5 


30 


21 


3 x 10 5 


1 x 10 5 


5.0 


0.45 


12 


7.1 


36 


28 


3 x 10 5 


3 x 10 5 


8.2 


0.86 


21 


12 


71 


50 


1 x 10 6 


1 x 10 5 


11 


0.98 


24 


13 


88 


59 


1 x 10 6 


3 x 10 5 


13 


1.7 


28 


17 


89 


60 


1 x 10 6 


1 x 10 6 


26 


2.9 


65 


32 


170 


110 


3 x 10 6 


3 x 10 5 


26 


1.1 


62 


32 


190 


96 


3 x 10 6 
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